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The Motion of a Solid in Infinite Liquid under 

no Forces, 

By A. G. Greenhill, F. K. S., Woolwich, England. 



The theory is sketched out in Thomson and Tait's " Natural Philosophy," and 
a complete solution for a solid of revolution is given in Kirchhoff's " Vorlesungen 
uber Mathematische Physik," IX; the treatment is developed at length in the 
"Motion of a Solid in a Liquid," by Dr. Thomas Craig (Van Nostrand, 1878), 
and in Halphen's "Fonctions elliptiques," t. II, chap. IV, and the subject in 
general has attracted the attention of a large number of writers, for which the 
Fortschritte der Mathematik may be consulted. 

The object of the present paper is to examine closely the elliptic function 
expression of all the dynamical quantities involved, and to explore the analytical 
field by working out completely the simplest Pseudo-Elliptic cases to serve as 
landmarks, utilizing for this purpose the analysis developed in the paper on 
" Pseudo-Elliptic Integrals and their Dynamical Applications," Proc. London 
Math. Society, vol. XXV, and carrying this out in continuation of the manner 
employed for a similar purpose in the papers on the " Dynamics of a Top " and 
on the "Associated Motion of a Top and of a Body under no Forces," Proc. 
London Math. Society, vols. XXVI, XXVII. 

The lectures delivered last year by Professor IQein at Princeton University 
have placed the analytical treatment of the motion of the top, and of Jacobi's 
two allied motions h la Poinsot, in a much clearer light and in a more elegant 
manner ; a similar treatment of the present problem will certainly prove equally 
valuable, but meanwhile the special cases, developed at length here, will serve 
as oases, so to speak, in the infinite region of the general elliptic function 
solution. 

Simple experimental illustrations of the motion can be observed in the 
evolutions of a plate or coin or bubble in water, or of a disc of paper or card- 
board in the air, as well as in the motion of a projectile or torpedo. 
1 
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1. The notation employed is that given in Basset's "Hydrodynamics," vol. 
I, Appendix III, and also in the " Applications of Elliptic Functions," p. 342 ; 
the body may for simplicity be taken as a smooth homogeneous solid of revolu- 
tion, having component linear arid angular velocities u, v, w and p, q, r; and 
now the total kinetic energy T of the body and of the surrounding infinite fric- 
tionless liquid stirred up by the motion of the body is given by an expression of 

the form 

T- \ P{u* + v*) + %R w * + $A(p* + <t) + i Gr\ (A) 

where P, R, A, G are constants depending on the shape of the body and on the 
density of the solid and liquid. 

The more general form of T assumed by Halphen, due to certain modifica- 
tions in the shape of the body, which still leads to Elliptic Function solutions, 
may be considered separately in its modification of the results, as also the effect 
of the circulation of the liquid, which may exist with a ring-shaped or perforated 
body. 

To realize practically the condition that no external forces act upon the 
body, even in a field of gravity, we may make the density of the body and of 
the liquid the same, so that the buoyancy and weight cancel, and the apparent 
weight is zero, as in the case of a fish and a submarine boat or torpedo ; and now 
the Hamiltonian equations of motion lead to 

p §^L — rPv + qRw = 0, (1) 

at 

p^ — pR w ^-rPu=0, (2) 

dt 

R <fo°_ q p u + pPv - 0> (3) 

A^- — rAq + qCr — wPv + vRw =0, (4) 

A^-—pCr + rAp — uRw + wPu=0, (5) 

Cut 

C^f- —qAp+pAq — vPu + uPv=0; (6) 

some of the equations being capable of obvious simplifications. 
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2. To make the memoir complete, it will be advisable to repeat to a certain 
extent the ordinary treatment; thus, multiplying (1) by Pu, (2) by Pv, (3) by 
Bw, and adding, 

^§»+p>i,*>+BSo $?■=*, (7) 

at at dt 

and integrating, 

P* (u? + v 2 ) + R*u? = F z , (B) 

where the constant -^represents the resultant linear momentum of the system. 
Similarly it can be shown that 

AuPp + AvPq + GwBr = G , (C) 

where G is a constant, representing the resultant angular momentum of the 
system. 

From equations (A) and (B), 

A (p* + q 2 ) = 2T— Cr* — Bw 2 —P (u % + tr) 



= 2T-GS-*~ + F*(±— -L)(*» - B?u?), 



and from equation (3), 

B*^ = P*(uq-vpr 

— P* (u 2 + v 2 )(p z + q 2 ) — P 2 (mp + vqf 

+ ( 2r _ a. _ £) r-w - ( g -«■*- )•, ( d) 

thus determining m> or i£w as an elliptic function of £. 

3. By the Principles of the Conservation of Energy and Momentum, T will 
remain constant, while F will represent a constant linear momentum in a fixed 
direction, Oz suppose; so that denoting by y lt y %1 y 3 the cosines of the angles 
between Oz and the axes OA, OB, OG fixed in the body, 

Pu=Fy u Pv = Fy 2 , Bw=zFy s ; (8) 
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and, introducing Buler's angles $, 4>, 4s 

y x = — sin 6 cos <p , <y 2 = sin 6 sin $ , y 3 = cos ; (9) 

p = sin <£ — sin 6 cos $ 4 , 
<2 = cos 4> + sin 6 sin $ 4, 

r— <£ + cos0 -Ji ( 10 ) 

so that 

P {up + wg) = i^sin ( — p cos <p + <? sin $) 

or <H __ (? — GrF cos 6 

~dl~ AF sin 2 

Split up into two partial fractions, 



= Fain* 6 ^ , (11) 

(12) 



c^ cfe efe 
where <H> X _ # — Cr.F 1 



(14) 



^ 2^.F 1 — cos d ' 

<&k _ g + O -F 1 

eft 2AF 1 + cos B ' 

and then ^. = r _ cos 0^ 

dt dt 

= 0-x)-f + f i - < 16 > 

4. Writing z for cos 0, or Fz for ifow in equation (D), then 

- =rf {. ( ,_ ir _(„.- ft ._^)^i-(*^)'}.(,., 



an* 



and according as the body is prolate or oblate, so is P > jR or P < i? , or an* is 
positive or negative. 

To distinguish these two cases, we take a = + 1 for prolate bodies and 
a = — 1 for oblate bodies ; so that 

F % 



TV 



= £(i~*). < 18 > 
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and now dz 

where 



dt =n^Z, (19) 



, = . ( ,_ 1) P_( M ._ tt ._£ ) * s i_(*»^y m 

(21) 



and 

.# ~ G—Crz 

dz AnF(l — z*)*/Z' 

Following the notation employed in the discussion of the motion of a top, it 
is convenient to put 

2T— Cr*—~ = An*D , (22) 

and g - PW _ „ mv 

so that „_ , . , W8 , n . /CrFz—G\* 



Z 



Z 



a(z>-l ){ z>-l-aD)-(V?^y, (24) 

= a(z*-l)(z»-l-aF)-(^=J^\ (25) 



5. Starting now with the general elliptic integral of the first kind, 

where the quartic 

^= az i + 45s 8 + 6cz 3 + 4dz + e, (27) 



and u denotes the elliptic argument (a different use of the letter u to that 
employed previously to denote a component velocity ; these two meanings will 
not however be found to clash) ; we make use of a theorem, supposed to be due 
to Weierstrass, but first employed by Biermann in his Dissertation, " Proble- 
mata qusedam mechanica functionum ellipticarum ope soluta," 1865, which 
asserts that we can put 

H^±-*)= F{Zl >^_ v z $ VZ > , (B) 

where 

F{z x , zg) = az\zi + 2bz& (z, + z % ) + e (a? + 4 2l s 2 + z\) + 2d (% + %) + e, ( 28) 

the second polar of Z; u x and % denoting the elliptic arguments corresponding 
to Zi and z 2 . 
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6. When the roots e a , e p , e y of the discriminating cubic 

4*-gf — g,= (29) 

are known, in which g % and g 3 are the quadrinvariant and cubinvariant of the 

quartic Z, given by 

g t = ae — 4hd + 3c 2 , (30) 

g, = ace + 2bcd — ad? — eb % — c 3 , (31) 

the above theorem (B) can be replaced by 

9 {th ± «*) — «a 

— g \ *S (zi—ZQ'Zi—ZqW (%—*?•%—*■*) =FV( g a— gQ-gji— g J^( g i— z 3- g i— z y) ] a ; ( 32 ) 
I 2(2i — z 2 ) ) ' 

where z , z a , Zp, z v denote the roots of the quartic Z=0, so that 

Z = a (z — z )(z — z a )(z - Zp)(z — a,) . (33) 

With the form of Z in (20), 

2 -r-3 a + z, + z v =0, (34) 

GV 
a (z + z a )(z p + z y ) + a (z z a + z„z y ) = — 2a — Z> — -^ , (35) 

a 2o Wv = a + jD — A z n i F z • ( 36 ) 

Thence, denoting by v x and w 2 the values of the elliptic argument correspond- 
ing to 

z x = + 1 and z z = — 1 , 

= sa{l-(z„ + z a )(z p + a,) + ¥« + Vv + z oWv } ± i — ^sw a — 

f (TV 2 ~ £ 2 CVjP 2 — (T 2 ) , ^ 

= §{a + 2a(z +z a ) 2 -2a-£- ^ + a+D- -^ ± ° ^^ }; (38) 
9{vi — v i )-e a = \a{z -\-z a f-^^ i ; (40) 
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* a ( Z + O 2 = 4A 2 n 2 F * + * fa + V ») — e « ' ( 41 ) 

(7V 

~~ ' +J?K — Vg)— e a . (42) 



4J 2 w 2 



7. It is convenient to employ the notation of Darboux in the corresponding 
motion of the top (Proc. London Math. Society, vol. XXVII), and to put 

G 2L Or 2B 



AnF M ' An M ' 
so that „ ,« , ws 1 r,, . fBz — L 



(43) 



Z=a(z*-l)(z*-l-aD)-4(™j^y (44) 

= aP-lW-l-aB)-4(**^\ (45) 

and §k. — 9 L— Bz , v 

dz ~~ ilf(l — z 2 )^* V ' 

Also s — s. ,-..,,. 

where ilf is a homogeneity factor to be determined in the sequel, and s is a new 

variable defined by 

u _ p ds 1 f* dz , . 

~M~J VS~™~WJ Vz y K > 

where S= 4s(s + a-) 2 — \{y + \)s — xy\*, (49) 

x and y being the quantities defined by Halphen in his " Fonctions elliptiques," 
t. I, p. 103; and when resolved into factors, we put 

8=4(s — s a )(s — s p )(s — s y ), or 4 (s — s^s — s 2 )(s — s 3 ). (50) 

8. Now putting 

Vl + v 3 =v, (51) 

(again a different signification of v to that employed at the outset) and denoting 
by o or s (v) the value of s or s (u) corresponding to u — v, equation (41) becomes 

$aM*(z + z a y = L* + o — s a 

= aNi, (52) 

suppose ; so that 

$M(z + z 1 ) = N 1 , 

IMiz. + z^-N^ 

$M(z + z 3 ) = N 3 . (53) 
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Taking Si>s 2 >s 3 , and employing relation (34), we find that when a == +1 , 

as for a prolate solid, 

Mz = —N 1 + N z + N 3 , 

Mz z = + Ni+N 2 —N 3 , 

M h = —N 1 -N t -N ai (54) 

while for a = — 1 , and an oblate solid, 

Mz 3 = + N 1 + N 2 -N 3 , 
Mz^ + ^ — Nz + N,, 
Mz 1 = — N 1 + N % + N 3 , 
Mz = -N 1 -N i — N s . (55) 

9. Now, rewriting the expressions for Z, in this new notation, 
Z= a (z — z )(z — a.)(a — z„)(z — z y ) 

= «(V-2 1 ^ j| ^^ 3 g 8 + 8 -^J-'a 

+ jg + ^ + ^ ~ 2iV ^ ~ gjg^ - 2NjNj \ . (56) 
so that, putting a= =fc 1 in (56) and (44), (45), 

H 1 2 — jf — + 8 ~i^~ + w ; 

( N? + N£+Nj 8 ^ja , Nt+....-2NjNj.... \ 

a \}-* — w 8 -j^- + w ; 

Subtracting, ffi-fiffi, __ 1Q BL^ 

M s M z 

or BLM=aN 1 N i N 8 , (59) 

so that » = » = fc *gi (60) 



with _^_ _ J2Z 

^J? if ' 



(61) 
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Multiplying equations (57) and (58) by a, and adding, 

-, o iV? + N$ + Nj N£ + N[ ±Ntz 2NJNJ - 2NjN[ - 2N?Nj 
M 2 "*" M* 

A & + D 

&M* - 2 [N* + NI + Nl — 2aZ 3 ) L Z M* 

+ (.A? + . . . . - 2N*Ni— ....)U+ AaNlNINI = 0, 
L 2 (M* — JV? — N$ — Ni+ 2aL 2 f = 4 (X 2 — aiV 2 )(Z 3 — aiVf )(Z 2 — a-fl?) 

= 4(* 1 — <r)(* 8 — <t)(« 8 — or) = — 2, (62) 

where 2 denotes the value of S in (50) when s = c ; and thus 

Jtf» = JV? + 2V? + Nl — 2aV — ^=1? (63) 

thus determining the homogeneity factor M. 
10. Equating the coefficients of z % , 

!k ' + D + i W = 2a + E + i W 

= 2 „:M±^ ; (64) 

so that p _ N? + iV. 2 + iV? — 2aZ 2 — ilf 2 

= 2a^=#; (65) 

^ i> = 2a-4=^4-4 Z2 -^. (66) 

But otherwise, equating the coefficients of a , 



Z 2 _B a 

ilf 2 ^ M 2 



«-t--0-4 ir?I = a + ^-4 ir .2 



iV? + — 2.AOT . . 

__ a i . if 



m 



(67) 



«+!>- 4Z 2 (iV 1 2 + iy|+iV3 a — 2aX 2 )— 4ZV— S + a(iy?+ • • - ■ ~^NjNj. . . .) (6g . 

if 4 ; 
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and 

N? + Nl + Nl = 3aZ a + Sapv , (69) 

Nf + 2iV|iV| = — 3Z 4 — 6L*M*pv + 9 JfV» — 2MYv ; (70) 

so that ,n_ (X 2 4-3Jf 2 H 2 — 4aiV(— 2)— 2ifV'« / 7 -i\ 

a -|- x/ — a j~ , ^ij 

an expression analogous to that in equation (214), p. 581, Proc. London Math. 
Society, vol. XXVII. 



We may also write 
and then by analogy, 



E= 




Mip'v 
2 L ' 


D = 




Mif'w 
2 B ■ 


w = 


Vi 


— V 2< 



(72) 



(73) 



where w = v t — v 2 , (74) 
again a different use of w. 

11. Another resolution of the quartic Z can be given, by means of the 

elliptic functions of the argument v 3 , which corresponds to the infinite value 
of z; namely, 

V(a)(s — ao)= ~ f% , (75) 

/v pu — pv 3 v ' 



Differentiating, 



and integrating, 



Va{z-zJ= ~ PV ° «"— - . (76) 

pu — pv 3 pv 3 —e a 

du (pu — pv s y 
= p(u — v 3 ) — p(u + v 3 ), (77) 

a/ (a) z = £ (u + v 3 ) — £(u — v 3 ) — X,2v 3 

— i P'(u — v 3 ) — p"2v 3 
* p(u — v s ) — p2v 3 



and squaring, 



t p' {u — v 3 ) + p' {u + v 3 ) , 
¥ P (u — v 3 ) — p (u + « 8 ) ' 



(78) 
az 2 = p (u — « 3 ) + p (w + v 3 ) + p2v 3 . (79) 
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1 2. Putting u x — u % = u in formula (E) leads to Hermite transformation 

p2u = — E-, (80) 

where H is the Hessian of the quartie Z, namely, 

H= (ae — ¥) z 4 + 2 (ad — 6c) z 3 + (ae + 2&d — 3c 2 ) z 3 -f- 2 (6e — cd) % + ce — d\ (8 1 ) 

Since b = in our form of the quartie Z (20), then corresponding to z — <x>, 

p2v 3 =-c, (82) 

or 6p2« 3 is the coefficient of — a 3 in Z; so that 

p2v s = ia + iD + i^ 



|a + |^ + 3 



if 2 



_ Nl+N j + Nj 

= ^+PK + ^)> (83) 

P («i + f>») = P2v 3 — 17 2 ; (84) 



and, by analogy, 



ilf 3 ' 



T>2 

P (*i — «0 = P2»« — -jut 5 (85) 



° r P(«i + ^) = |a + l^-|-^, (86) 

i? 3 

Pfa — v 2 ) = i« + iD— -^j (87) 

while from (72) and (73), 

W fo ~v z )= *j> (90) 

So also Hermite's formula 

P'2w = ^-, (91) 
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where Q denotes the sextic covariant of Z, leads for the infinite value of z to 

?L__s = coeflScient of 4z in Z, 

Vf*H= Wa™. (92) 

We notice that v 8 is a fraction of the real or the imaginary period, according 
as a is +1 or — 1 , or as the body is prolate or oblate ; and, sometimes, by 
retaining a, we are able to state the results in a general form, suitable for all 
cases. 

13. The annexed diagrams are intended to illustrate the various cases which 
may arise, and to exhibit to the eye the separation of the roots of the quartic Z ; 
the curves on the left hand are the typical graphs of the function Z, while the 
concentric circles on the right show the correspondence of z on the outer circle 
with the elliptic argument u on the inner circle, the shaded portion representing 
the limits of the actual variation of z. 

Case I applies to the prolate body for which a = -f- 1 ; and all four roots of 
Z are real, and arranged in the order 

00 > Z! o>l> 2 3>2>z 3 > — 1 > esj > — <», 

so that %, o 3 denoting the real and imaginary periods of %>u, and the letter/ 
being employed to denote generically any fraction, 

*i =/<»3> «a = °i +/%> v s =/%, 
v or w = v 1 ±v 2 = a 1 -\- fa 3 , 
u = a 3 +/a 1 . 

Case II, an oblate body, a = — 1 ; all four roots of Z real and arranged in 
the order 

00 > 1 > »3 > a > 2 2 > % > z > a > — 1 > — °°. 

*>i, Vz, v 3 , v, w = /a s , 
w = 6) 3 +/% or /%. 

Case III, an oblate body, a = — 1 ; two real roots of Z and two imaginary : 

°°>1>Z0>3>%>— 1> — », 

v lt v it v„ v, w= f<4, 





— 1 



14 Grbenhill : Motion of a Solid in Liquid. 

Cases IV" and V represent the conditions for an oblate body when the 
kinetic energy assumes the more general form discussed by Halphen and others, 
in consequence of which Z takes the more general form, with a = — 1 , 

Z=-(z*-l)(z* — 4 ^-l+D)-4(^=^) 3 



= - (z 3 - l)(i»-4ta- 1 + B) - 4 (^y^J; (93) 

and now in 

Case IV, QO >0o>3i> 1 >2 2 >3> 25 3>— !> — °°> 

v 1 = a 1 +fa 8 , v % and v 3 = fa 3 , 
v and iv = g> x + /% • 

Case- V. 00 > 1 > z 3 > a > z a > — 1 > z 1 > a >— 00 , 

«! and «? s = /w 3 , v 2 = u 1 +/o 8 . 
-w and to = Wi + /fo 3 . 

The term & may arise from the external shape or from circulation in the 
liquid when the body is perforated or ring-shaped (Basset, Hydrodynamics, 
I, §193), but the absence of b simplifies considerably the elliptic-function expres- 
sions. 

14. Putting z= ± 1 in (77) and (79), 

W (-a)^ = J 0^ = ^,- Vt )-*( V ,+v,y, (95) 

a = fp(v 1 — v 3 ) +p(v 1 + v 3 ) + p2v 3 , (96) 

a — ^{v % — v 3 ) + fp(v 2 + v 3 ) + $>2v 3 . (97) 

ThenCG 9 K - O = * a - i *>2*, 3 + V(- «) ^^ , (98) 

7? — 7 

$>K + = ifl — it>2P 8 — V(~ «) ^ » (") 

p I r 

pfa — «*) = ia— ■ip^j + V'C— a)— jg— 1 ( 10 °) 

P I 7" 

p(v % + v 3 ) = ia — ip2v 3 — */(— a) M . (101) 



and 
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Again, putting z= ± 1 in (78), 

P' fa — v 3 ) — p'2v 3 ~ Wa \fpfa — v 8 ) — p2v 3 \, (102) 

&'fa — v 3 )= 2 N /a|-^g +ia~fg>2v 3 + V(— a) ~ I 

= W«^-iVa^-2Va^ + 2i^^, (103) 

£>' (»i — v s ) + jp' (»! + v 3 ) = 2V« {P (»i — v 3 ) ~ <p fa + t> 8 ) } 

= 4 *T' ( 104) 

P , fa + v 3 )=-Wa^ + ltVaE + 2Va±r i + 2i^ w ^. (105) 



'(* 8 --^)= 2Va^ + iVc^+2Va^-2t^±^, (106) 

/(^ + ^) = -Va^-W^-2Vfl^-2i^. (107) 



Similarly 

M % M 2 M 

Thus 

9 fa — «») — P (v a + B,)= 2V(— «)-^, (108) 

d r i> 

&>' (»i — «*) — *>' («» + » 8 ) = Wa ^gf + 4 * Y ; ( 109 ^ 

so that 

, fa - v 3 ) - ,>» + vj ( L 

1pfa — v s ) — pfa + v 3 ) M v / 

and squaring, 

SP(«i + ^)+P('«i — v B ) + pfa + v B )=— -jp — 2*/{—d)-yj:+a. (Ill) 

P(«i — v 3 ) + pfa + v 3 ) = a — &2v 3 — 2*/(~- a)-gr ; (112) 

so that , . , _ I? /, - „\ 

fK + » 3 ) = F2v 3 ~--^, (113) 

agreeing with equation (84), 
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Similarly, 

P(»i + v s) — V (?* + *») = W(— «)-g=, 

P' (vi + v B ) + p' (v 2 + »,) = — 4 V« ]gr — 4 * y 

so that, squaring, 

jB 2 B 

P{vi — v t ) +9(vi + v 9 ) + p(v a +v 8 ) = — -jp — W{— a) -^ + a. (115) 

But P(»i + ^)+ P(»,+ «,) = fl-A-2V(-o)y 1 (116) 



f»(«i — «i)=l»2t>8- -jgr. ( 117 ) 



so that .. , , B 

agreeing with equation (85). 

15. Again, from (110), 

^-P'^-^V^-P'^-Vs) jL (U8) 

Pi^ + v,) — 9 («! — « 8 ) if^ V V ; 

and 

tf> (*i + O — P K — v s) = ^2« 3 — ^p- — £a + i f1v z — s/(— a) —jg— 

= lE-^{-a)E=Lk, (119) 

therefore 

V f fa+*%)= 2 (* jf ~ *Sa)\P( v i + V2) — P( v t — v s)\—P'( v 2— v 3) 

=(4- va ){^- v( - fl ^l 

= **T« (120) 

as before in (89). 

So also, from (114), 

y> fa — v> ) — ff fa -j- p 8 ) _ 1 — l/fa + q,) — l/fa + pg) 
f(«i — «») — !? (»8 + « 8 ) S P Oi + *s) — P (»8 + *>») 

= -i^+Va, (121) 
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and 

732 P i r 

F(»i — « b ) — pfa+ t> 3 ) = jf>2» 3 — -gy — i a + | p2w 3 + V(— a) ^ 

= 4Z?+V(-a)^±^; (122) 

therefore 

agreeing with (90). 

These verifications are useful in fixing the signs of the various expressions, 
in the ambiguous cases, of frequent occurrence in these calculations. 

16. Returning to (75) and putting in it ti = v l , and z= 1, 

V«(l-2 )= ~ f ' V * , (124) 

pv 1 — pv 3 

SOthat Va(l-,) = , -^3(^-^1) (125) 



and similarly, 



(f»»i — P» 8 ) (#> M " tw 8 ) ' 

^ a (-l- 8)=a -t^fr"-*) (126) 

(#>»8 — pt'sHjw— pt? 3 ) ' 



Also, from (24), (25), (44), (45) and (77), 

— GrF _L — B_ 1 f/n f/r, 



(127) 
(128) 



2AnF M 2V(— «) (PV— F^) 2 ' 

ff + CrF _ L + B _ 1 g>'?; a ^ 

24^ if 2V(~ a) (JW 8 — P^s)* ' 

so that in (14) 

d4£_ . a — GrF 1 
di* * 2-Aw.F 1 — 3 

_ 1 P'vi (pw — fw 8 ) 
(p»i — F»») (fw — pvj 

§Wi — fPV 8 PU — JWx 

= K(«i + «s) + i £(»i — «i) — fo 

- K(« -th) + K(" + »i) + fri- (129) 

3 
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Similarly 

d^i _. i G+ GrF 1 



du 2AnF 1 + z 

— 1 P' V 2 I 1 V' V Z 

•f 'S" 



pv 2 — pv 3 fu — fv % 
— KK + i's) +i£(% — «s) — ?«» 
-K(« — ^ +*£(« + «*)+&• ( 13 °) 

Integrating (129) and (130), 

4' 1 i=H^v 1 + v 3 )^Uv 1 -v 3 )}nt+ilog^±^ y (131) 

lW=*K(^ + ^)+f(«b-^)f«« + ilog|^±^; (132) 

and adding, 

4. = i Qnt + * log ^^ , (133) 

where <?==£(«i + «») + £ fa — «b) +£(«* + »s) + £(»i — «s). (134) 

By a theorem of Elliptic Functions, 

pu + Slgfc +5^ = ilogffi + ^) + ft, (135) 

°<5(w — «j) 6 (m — v%) 6(u — v) v ' 

where Vi + v^ — v, (51) 

and . _ tan _i (ggi — yea) j/g 

*V v a (pu — p»i) — *>S (pw — &vz) 

= sin_1 -VT7 w~^ v n • ( 136 ) 

so that we can put 

ty = $Qnt + Zi + i logi^±^ , (137) 

and we have now practically added the parameters v t and v 2 of two Elliptic 
Integrals of the third kind into a single parameter v. 
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17. We now introduce the Elliptic Integral of the third kind, with elliptic 
parameter v, in the standard form we shall employ and denote by I{v), namely, 



j( w ) = ^^)(fL~5)ryzzS 



(s— a)VS 
(u — v ) 
6 (u 4- v) 



it- log pp=*l £•*-'£)*<«•*, (138) 



with 

/ds u nt /ioo\ 

»/S MM v ' 

where P(v) or P is a certain function of v, to be defined hereafter, so chosen as 
to cancel the secular term when the integral I(v) is pseudo-elliptic, in conse- 
quence of the parameter v being an aliquot part of a period. 
Now from (138), 

i log *( u + v ) = U(v) - n%v -H-£- nt, (140) 

6 [u — v) M 

and thus 
fy = iU («i + «j) + ? («* — «»)■— £ (»i + »») 

+ ?(«i. + «s) + f (»i -»s) -?(»i + «*) - »'-^-| nt + il(v) + £*. (141) 
But, from (99) and (100), 
£ 0i + «*) + £ (v$ — v s) — £ («i + »*) 

= V ]p(«i +« s) + P fas — «s) + P (*1 + »») f 

+ I a — lp2v 3 + V(— «) + + P2v 3 — -jp \ 
= v{a+2V(-a)A_^| =4 + ^ a - (142 ) 

and similarly, 

£ («*+ vs) + £ K— «s) — £ K + *>a) = • ^ 

To settle the ambiguity of sign, we may also employ the formulas 



£ (vt+ vs) + Z(vi — v 3 ) — ?(«i + '» 2 ) = i 1 rf — Va. (143) 



y /• . \ . y / \ >• / 1 \ — 1 P' (^1 + «s) — P' (*>g — ^3) 

g fo + , 3 ) + g K - ^ - g {V1 +„) = - ^ ft + *) - y fa- *\ 



(144) 
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and thus 

L ,, P 



^ = (* ~M ~ * * ~m) nt + U ^ + **' 



and, dropping the factor »', 



or 



* = ^=jpk* n< + /(«) + £, (145) 

4,-pt = I(v)+£, (146) 



where ^ L — j P(v) 

~n M 



(147) 



a different^? to the component angular velocity about OA in §1. Without this 
preliminary determination of pt, the secular term in ■vj', the consideration of 
the pseudo-elliptic solutions would be hopeless. 

18. The next chief object of the present paper is the discussion of the 
Pseudo-Elliptic cases which arise when the parameter v is made an aliquot part, 
one n th , of a period, so as to be of the form 

v = 6)l + ^3, or B2S, (148) 

n n 

according as the body is prolate or oblate ; and when P (v) is at the same time 
so chosen as to make I(v) an inverse circular function ; the preliminary analysis 
will be found in the paper on "Pseudo-Elliptic Integrals," Proc. London Math. 
Society, vol. XXV, from which the results required in the sequel will be taken. 
In such cases it will be found that we must take 

(149) 

JXL tO 

and .•r/.. l _n_<3(« + «) J i. (150) 



where 



%1 


M 


)- 


£v- 


5 


:i(v) 


= : 


Hog -A- 
6 [u 


+ «0 e *« 

— v) 


V 


= 


n 


, s: 


1y\ 
~~n~' 



(151) 



and now il(v) is the logarithm of a function, analogous to the sn, en or dn function, 

which is considered in Halphen's " Ponctions elliptiques," I, p. 224, the function 

being the n th root of a rational function of $>u and p'u. 

Further, by taking 

L = lP(v), (152) 
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the secular term pt is cancelled in the equations and all trace of elliptic trans- 
cendentalism is eliminated, so that the various curves described by points in the 
axis of the body are, relatively to the moving origin 0, purely algebraical 
curves ; these special cases are interesting to discuss and to represent diagram- 
matically and stereoscopically, like the analogous algebraical Spherical Catena- 
ries and G-yrostat Curves, investigated in the Proc. London Math. Society, vol. 
XXVII, and drawn by Mr. T. I. Dewar. 

19. The curve described by the projection of the moving origin on a 
plane perpendicular to Oz is intimately associated with the cone described by 
the axis of figure 00 round Oz; for, denoting the coordinates of the projection 
of by a, /?, and the advance of parallel to Oz by y, then, according to the 
equations given in Kirchhoff s " Vorlesungen," p. 240, 

Fa= frVL + frML + p,™, (153) 

F P=- a '-dp— a >-dq-- a *-W' (1M) 

w dy dT , dT , dT ,,__, 

F jT~ u -^ r- v -g— + w -*— ; (155) 

at da dv dw v 

where u, v, w again denote the component velocities of in the directions 
OA, OB, 00) and a lf a % , a g denote the cosines of the angles between O'a 
and OA, OB, 00; and (3 lt /? 8 , /? 3 the cosines of the angles between O'fi and 
OA, OB, 00; Oa, O'fi, Oy being three fixed axes in space, drawn through a 
fixed origin O, O'y being drawn parallel to Oz . 
Expressed by means of Euler's angles 6, <?>, 4, 

a x = cos 6 cos cos 4 — sin <£> sin 4, 
a 2 = — cos 6 sin <£ cos 4" — cos <£> sin 4', 
a 3 = sin 6 cos ■fy; (156) 

^1 = cos 6 cos $ sin 4 1 + sin <£> cos 4 . 
@ 2 = — cos 6 sin ^> sin 4 + cos <£> cos 4 . 
/? 3 = sin 6 sin 4; (157) 

* " 

while p= sin<£>0 — sin 6 cos <p4/, 

g- = cos 4>0 + sin 6 sin $4> 
r = $ + cos $4>; (158) 
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so that, after reduction, 

Fa = A cos 1^0 + ( Gr — A cos 64) sin 6 sin ^ . 

^/? = J. sin 4>0 — {Gr — A cos 04>) sin cos 4<, (159) 

j 
F ~j— = a An 2 sin cos cos 4 1 1 

F-¥- = aJ.n 2 sin cos sin 4 1 . 

jF *L - *1 + oj1w» cos 2 0. (160) 

20. Changing to polar coordinates p , m in the plane O f a(3 , so that 

a = p cos or, /3 = psinnr, (161) 

then 

>(a + pi) = Fpe m — { Ad — i sin ( Gr — 4 cos fr£) \ «**, (162) 

- d J|,oo.(*— ) = ^=-^^, (163) 



sin0 ' 

ify sin (4- — nr) = ( Cr — J. cos 0^) sm 
GrF — G cos 



and, dividing, 



so that, from equation (146), 



. _ x GrF— G cos 
*9 
and nr and 4* depend upon the same elliptic integral I{v) 



(164) 



F sin 
Thus, squaring and adding, 

* 9 — F*{l — z*) 

— A*n*a {aE+ 1 — z 2 ); (165) 



•»<+— >=-toz ! ' (166) 

- = ^ + tan ^iV^ 

, . • _i GrF— GcosB 
= 4* + sm * =* — = — 5 — 

.F a p sm 

^^"TMpfcl 1 (m) 



w=-pt + I{v) + £ + sm — g . n e , (168) 
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We find in fact, by differentiation of (167), 

dm _ d^ , d_ . _ x GrF— Gz 

dz ~~ ~dz~ dz Sm AnF*/\a{l— z 2 )(aE + l — z*)\ 
_ g —GrFz J_ 
AnFil—z*) */Z 

Or F— Gz / — G _z_ z \ 

AnF^{a{l — z*){aE + l— z % )\ \CrF—Gz _ 1 — z 2 + aE + 1— eV 

i " V^ 

V^a(l— 2 a )(aJ^+l— z 2 )} 

f G — GrFz G , GrFz — Gz 3 , GrFz — Gz 2 



+ 



GrFz — Gz 2 ] _1_ 
"•" AnFiaE + 1 — s 3 ) J -^ 



. 4nJ'(l — z 2 ) ^Ini? 4rci? (1 — a 2 ) ^ AnF(aE + 1 — z 3 ) j V^ 
GrFz—Gz % 1 



AnF(aE + 1 — z 2 ) V^ 

~ 2 if (ai? + 1 - z 2 ) ^Z* (169) 

Now if ic l9 w 2 denote the values of the elliptic argument u corresponding to 

aE+1 — z 2 = 0, (170) 

it follows by Abel's theorem and the theory of elliptic functions that 

«>1 -f Wg= Vi + « 2 = v. ( 17 1) 

So also if t x , t % denote the values of u corresponding to 

aD + 1— z 2 = 0, (172) 

then t y — t 2 = v 1 — v 2 = w. (173) 

21. Prom equations (125) and (126), 

sin 2 = 1 — z 8 = — a <A 3 (fm — ggiXgM — pv») 

(9®! — «»»,)(f>cjj — pv t )(pu — pv 3 f 

— /7 g (M — «i)<3(*« + «i)6(m — p»)6(m + «i») / l74 x 

<S 2 (w-« 3 )<5 2 (w + « 3 ) ' lW4j 

*» = g , (w + gi| g $ M + *> e**; (175) 

6 (m — vj) 6 (u — *? 2 ) v ' 



and from (133), 



so that by multiplication, 



sin W = y g gi?i± gi) ^ (" ± *■) e* -. (176) 

6 («* — «? 3 ) 6 (w + ^3) 
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Or, in Klein's manner, with the stereographic projection, 

tan j Q* 1 = C> PE± *> g ^ ± *j 6- «*. (177) 

6 (it — « x ) 6 (u — v % ) 

Similarly we find 

wi —(jn g ( u ± w l) gJM + «>») g*fl* ( 178 ) 

so that the expressions in Halphen's equations (52), F. E. II, p. 164, appear to 

require correction ; however, the curve of (a, (3) is intimately connected with the 

cone of (0, 4')> and the two are pseudo-elliptic, and even algebraical, at the same 

time. 

As for u, v and p,q, they depend upon the same elliptic functions as 

Euler's angle $ ; for 

P ( u + vi) = — F sin e~*\ (179) 

p + qi = (— sin 0j/ + id) g~** 



= — n t— nx e" 



-04 



(180) 



V(l — a*) 
which are pseudo-elliptic together with <£> ; and this is the case when 

V\ — ^2 = h — h = w ( 18 1) 

is an aliquot part of a period ; but these cases are not so interesting from the 
dynamical point of view. 

22. In the pseudo-elliptic case of the motion of a prolate solid, when the 
parameter 

V = tal + **L, (182) 

ft 

where n is an odd integer, the expressions for 4 1 and w must be of the form 

+ — P=— cos * (l — «?)*» e-Vfe — »•* — »,) 

= ^-sm- 1 (1 — g 2 )^ V(% — a.« — «.), (183) 

and 

« ~ J* = 1 ""- 1 Jrf> " X + (f +" 1 -'^" + — ^-«-»-^ 

= T sin- ^ gK " 1+ ( ^+;-+^ </(«,-..« -%)• (184) 
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But with a parameter 

w = G)l + _^L, (185) 

where n is an integer, the Zi and z % or z 3 must change places, as well as the 

cos -1 and sin -1 , taking -=- as positive at the start. 

With an oblate body the parameter v is a fraction of the imaginary period, 
and when 



v=-^-, (186) 

n 



where n is an integer, we must have 

^—pt—— cos (! _ ~#fi WCg-gi.z-gg) 

= A- sin- ^- ± %'_+ yV ■ + jr - ^-...-^. d«) 

with a similar expression for m — pt. 

When, in the motion of an oblate body, the parameter 

«=r— 3 , (188) 

n 

where n is an odd integer, the quartic Z will not be divided into quadratic fac- 
tors, and we must have 



4> — pt = — cos J ! \- U^ n 

* n (1 — z 2 )* n 



(i — zy 

= ± rin-i^-'+^r'+'-'+^-V^ (189) 

n (1 — 2 3 ) sn v 

with a similar expression for w — pt. 

23. Having now chosen a simple numerical value of n, for a case worked 
out in the paper on "Pseudo-Elliptic Integrals," and having expressed the 
s l , &,, s s in (50), or Halphen's x and y in (49), in terms of a single parameter c, 
as well as cr, -</ — 2 and -P(«), and having chosen an arbitrary value of L, the 
quickest practical procedure for the determination of the coefficients H, K and 
4 
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I, J, appears to be effected from the identical equations obtained by squaring 
and adding, and also from the differentiation of the two equivalent forms for 
4 and vs. 

Thus in consequence of the relations 

-| = WS, (19) 

c^_ n G- GrFz . . 

we have dfy „ L — Bz 1 , . fi v 

W~ ~M\\ — #) VZ' [ } 

and putting 

+ —pt = z, (190) 

d% — 2 L ~ Bz _J JP L 

dz M{\ — a a ) \/Z n »/Z' 

(l-z*)VZ^=2 L - Bz -^(l-z») 

v ' <fe M n v ' 

if 2f ^ Jf v ; 



So also 



(aff +1 — *)*/Z d ( a -rt 

24. Begin with the simplest pseudo-elliptic case of all, in which 

v = a a and tfv=0; (193) 

so that, from (72), L — 0, (194) 

or E—0. (195) 

When X = or 6? = , 

~dY = ~ 2 w (i-z*)vz' (196) 

d& _ 9 ^S_ _3 nqrl 

<fc "" M (aE+l— z*)VZ' Kl > 
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where Z is a quadratic in a 2 , so that taking s 2 as independent variable, ^ and nr 
are non-elliptic, and integrating, 

^ = * 8in " 1 Af^)' (198) 

- = i "°" APTT37 )' (I99) 

where A denotes the discriminant of -2T. 
For a prolate body, a = + 1 , 

Z=(fi-*){4-#), (200) 

. _, /(#+l — z\){z\ — z-) 

- eo^-i./ W-l-iQW-^) • (202) 
which may be written 

V(4 — 4) sin &** = V(l - 4V(4 — * 2 ) + *V(zf - 1)V (4 ~ A ( 203 ) 

V(4 - 4) ^ ^ roi =V(4 — i— ^ V(4 - * 2 ) + »V(^+i— 4V(4 - * 2 ), (204) 

and these equations prove that the cone described by the axis 00 round the 
fixed direction Oz is a quadric cone, while the curve of (a, (3) is a conic section, 
an ellipse. 

For an oblate body, a = — 1 ; and in Case II, 

Z = (4 — s?)(s? —. »4) (205) 



sin * 



7 (1 - 4)(4 - **) (206 ) 

V( Z §_ z f)(i-^)' ^ Ub ' 
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a -i J(4-L+l)E-i) 



*- C0Br v (4-^(^-1 + -^) 



B in-» ,/ W-l + ^)W-^) • (207) 



or 



V(zl — 4) sin ^=V(1 - 4V(z 2 — *S) + »V(1 - z!)V(a§ - 2 s ), (208) 

V(^-20^P eWi =V(^-l + ^V(z a — 2?) +*V(zf-l + ^V(4 — a?), (209) 

so that OG describes a quadric cone, and (a, fi) describes an ellipse, as before. 
In case III, with L = , we must put 

Z=(z\ + ^)(zl-^), (210) 

and now 

^(z\ + zf) sin W =V(1 — aJV(ef + z 2 ) + »V(l + 4W{4 — a 3 ), (211) 

V(z\ + 2l)~ p e ^=v/( z l— i + ^)v( 2 I+s 2 )+^(- Z f— 1 + #)V(4 — a 3 )- (212) 

The curve of (a, /?) is still a conic section, and OG describes a quadric 
cone relatively to Oz, but the azimuth ^ oscillates between 

± cos-Vi=4^= ± sin-V^44^. ( 213 ) 

z l + 4 4 + 4 

obtained by putting s = . 

Next, with E = , we have 



^: 



= a(s ._l)i_4(^L=^)". (214) 

With an oblate body, a = — 1 , and ^ is always negative, so that no solu- 
tion exists. 

But with a prolate body and a = + 1 , 
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. , ._ y( i - 2 i + 2 ff- g 
y( i+2 4- 2 f- o 

_ cos _____^ 

= i s in-i-^, (216) 

1 — s 

where p L A dz , „ /01 „ n 

» = ir .^- s - = »VZ, (216a) 

as may be verified by differentiation. 

Also, with E= 0, 

dm _ Bz — M 1 , 7 . 

"a¥ _2 F(r=^) V^' ( } 

so that A* , 4_ 9 £ 1 _ 9 „ * / 91R \ 

m + ,J,= 2/rf + -|-; (219) 

and thus , _ 1 V( 1 — 2 lT + 2 ^ _z2 ) 

t» »£ = COS 1 77T-7T5 9\ 

r */2*/(l — O 

Q1T1 -I 



= sm" 



V2V(1 — 2 2 ) 



= * sin ~ 1 ^' ( 22 °) 

25. With r = or i? = 0, the case of no rotation round OG, 

w = a a , (221) 

and Z= a (rf - !)(*- 1 - aD) - 4 ^ 



= a («» - l)(* — 1 - atf ) - 4 I* z z , (222) 



if 2 

J? 

M 2 
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and now equations (19) and (12) can be integrated by means of the ordinary 
Jacobian elliptic functions of the second stage. 

When i? = , then either N± or N 3 = , according as the body is prolate 
or oblate ; so that 

i = V(s 1 — a), or — */($ 3 — o), 

and M = >/(<? — * 3 ) — V (<* — <%) . or V(*i — c) + V(s z — a) ; 

1 4-aF- \ y(<>— ** ) + */(« — *9) W r i y( s i— <*) — */(*» — a ) V. 

i -ran- ■j V ( (T _ S3 )_ V ( <y _ Sg )} ° W(*i + <r) + V(*,-<r)J ' 
and 

# = |V _ f v(%— ^) — V(^i — gg)) 8 ~ir s g _ W(*i — s s) + vxgi — gg) | 'i 

L W(ff — * 8 ) — /s/(cr — s 3 )j JL W(<T — S 3 ) — \/(<T — s 2 )J J 

or __ |>_ f V( g l — g :0 — Vfe — gg) r~ir^ — W(*i — g s) + V(a g — gg) )'-) 

L W(*i — <r) + V(«, — o)) JL (V(*i — <0+v(*i — *)J J' 

With an oblate body we shall find that this makes .5 = 0, 5'= 1, _Ki=0, ... ., 
while in the case of a prolate body the parameter must be of the form (182) for 
similar reductions to take place ; and for the parameter in (185) we shall have 
H= iT=i\/2, when B andiV r 1 = 0. 

26. Passing on to the next simplest pseudo-elliptic case of bisection of 
a period, by taking a parameter 

v = a 1 +hG> s (223) 

for a prolate body, and writing % for 4- — pt throughout, the result must be of 
the form 

£ = 2Cos - ^ W {z — z.z 3 — z) 

= * sin_1 K l-^ v( - z ~ z% ' z ~ zl) (224) 

with K> + H*=l. (225) 
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The solution must be built up from the associated pseudo-elliptic integral 

TV _1_ 1 N i / , (c+C 8 -«)-2(c + C 2 ), 



c + c 2 * 

= i COS" 1 — ^ 
c -4- c 
in which 



= J cos" 1 , V * ; (226) 

C + <f 8 



Sl = (l + cf, s % = c\ s 3 =0; (227) 

P(v) = 1 , ff = c + c a , V— 2 = 2 (c + c 2 ). (228) 

Then, with the preceding notation of (54), 

Nl = II— 1— c, 
Ni = L> + + c, 

Ni = Z a + + c + c 2 , (229) 

so that 

Ni + N?Ni = L\ (230) 

.tf? - Z 2 = c + c 3 , (231) 

2% + Nj =21?— 1; (232) 

and from (59) and (63) 



if ZJf a ' 

= Z 2 -l + c + c 2 - 2c+2c2 



(233) 



Z ' 

iJf 2 = 2Z 2 — Z + (Z - 2)^?, 

Mz % =z+N l +N % — N 3 , 

Mi = — N 1 — N z — N s , (235) 

arranged in the order 

°°> zo>i>%>2>22>~i>»i> — °°; 
with N.KN.KNs. 

To determine the coefficients Zf, ZT 1( K, K x in equations (224), differentiate 
with respect to z ; then since 

(z — Z){Z S — z)=2S 2 — 2-^-3 + 3 2 3 , 

iV 
(z — z 2 )(z — z a ) = s 2 + 2 -^ z + zfa, (236) 
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d % — 
dz 



■'WVC-4-+-) 



2 



E , * IF , 2s 



+ s^ + 



T 2 + Zo% 



/ 4- A 

/f ( z2+2 f *+*«.)(!-*) + («+ §)(^+ir i )(l-^)+22(^+/f 1 )(^+ 2^ Z + z 1% ) 
~ * -■-— — (Hz + ff^l — z^Z ,(237) 

each of which can be equated to (191), so that 

= #(s 2 -2§z+%)(s*-l) +(2-§)(^ + J ff 1 )(^-l)-2(^ 3 + i^)(V-2 Js + %) 

= 0. 2 4 + (-^ + f ^^{-(H^ff+sl^^H (238a) 

(m + H 1 )[— sr *-i- B * + ^^) 

= -K(*+2*±z + %z 3 )(z 2 -l) -(z+ ^){Kz + Eft*- 1) + 2 (K* + E,)(f + 2 1? * + *,%) 

= O.Ji' + ^.Ki — ^i)ss» + j(2 + 8 1 a,)Jr+8^ J Ki} S ?+....; (2386) 

and equating the coefficients in these identities gives us equations enough, and to spare, to determine 
-ff. -Si, K, K x . 
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Putting z = =fc 1 , gives the relations 

2 (fl ~ ^ ^i~ = _ (*- K K l - 2 jfr + % ) • (239) 



which are useful forms to serve for verification. 

"Working in this manner we find after considerable reduction that 
H _ 2? — iV^ if _ Z 2 + ^i^ 2 





JT 


N 3 ' H N z 


and, squaring 


and adding 


in (224), 


so that 




1 K , H _2L* 

EK~ H ^ K N 3 ' 


or 




2HK=^, 


and 


K 
H 


~:r- 2 -ivr' 



so that 



Also iT, _ (Z-lf-iV^ 

*i - -(z-i)*--#ir, 



(240) 
(241) 



(242) 



K*—H*=z^; (243) 

If 

2ir~ ' 

m— 1 ?— N i N z (244) 

^ ~" 2i 2 V ' 



(245) 



so that &! , Zi 9 JVi^g ^246) 
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as is verified by squaring and adding ; and all other relations are also found to 
verify. 

27. In the case of the oblate body with parameter 

« = *©,, (247) 

we must take 

i — 1 Sz -4- S-i , i \ 

% — i cos — ^ */(z — z 1 .z — e,) 

= I sin" 1 K [^_g x V(z 3 — z.z- z ) , (248) 

and build up upon the pseudo-elliptic integral 

I(i a \-i Al + 2c)(e ± o + c 2 ) — 2 (1 + 2c)(c + c 2 ) , 
1 ( 4 o 8 ) _ ij j- - - - ^ s ds 

= 3 sm 1 -5U v , / g [ 

1 _i(l + 2c)Vs /o^n\ 

= icos 1 ^— p- — PV, (249) 

s + c + c 4 

so that, with the same s lt s z , s s , we take 

P(v) or P = l + 2c, <? = — c — c 2 , V— 2 = 2(1 + 2c)(c + c 2 ); (250) 

and now, with (52), 

Nf = 1 + 3c + 2c 2 — X 2 , 
JV| = c + 2c 2 — Z 2 , 

2v? == c + c* —I?, (251) 

so that N* + N£= p 2 — 2£ 2 , (25 2) 

2V| + X 2 = c + c 3 , (253) 

JV?^f — i 4 = P 2 iV| ; (254) 

and from (59) and (63) , 



a ' 



M LM 

jf 2 =l + 5c4-5c 2 -Z 2 - 2 ( 1 + y + * 2 ) 



(255) 



= p^ + NI— 2~{Ni + U). (256) 
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Proceeding as before, we now find 

H_ „ N X N S + L % JT _ N^—L 2 
K ' PN S ' H ' PiVg 



(257) 



and H % — IP = 1 , (258) 

2HK=™*, (259) 

H 2 + K 2 =^, (260) 

2Z 5 ' 
^ = ^ ~ £2 , (261) 

s, _(L ~ Pf + jvijy 8 
if ~ pjf 

# pjf 

and now all the conditions are found to be satisfied. 



(262) 



28. Both results for the prolate and oblate body with parameter 

» = Ui + io 3) or io 3 , (263) 

can be incorporated in the form 

Z = i cos -1 *_g 1 V(z — z . z — z 3 ) or V(z — % • z — %) 

= J sin- 1 Kz + ^ s/(z — z x .z — z») or >/(z 3 -- z . z — z ) , (264) 

1 — z* 

with J EP + a J ^ 2 = l, (265) 

H> — aIP=—aW*, (266) 

2HK=^1, (267) 

If we try to cancel the secular term pt by putting 

I = l?=i or i(l + 2c), (268) 

we find N?= — I — e, or Ni= — $, (269) 

thus leading to imaginary results. 
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As verifications we may take the numerical case for a parameter v = a 1 -\-la s , 
worked out on p. 534 of "Les fonctions elliptiques et leur applications," Green- 
hill-Griess, 1895. 

For a parameter v = ia 3 , take, as a special case, 

JB=0, N 3 = 0, JVi=l + e, N 2 = c, 
Z=— V(c + c 3 ), M— V(l+ 2c){v/(l + c) + Ve\, 
P- = L — \P - _ V(l + c) + \/c . 
n ilf V(l + 2c) ' 

and now 

♦-^'"-^[■•- Mi yi 

y/(l+ C ) + VC 

-^sia- 1 _^lL±|^_V[(l + 2c)|V(l+c)— VcP— z 2 ]. (269a) 

29. Proceeding to the next case where the parameter 

« = Oi + i<3 3 » or io 3 , (270) 

both results for the prolate and the oblate body can be incorporated in the form 

^ = ^-cos - _J ' — W(z 3 — z.z — Zj) or »/{z — z 2 . 3 — z x ) 

= i sin ' _ l J i — - V(z — ^ • z — %) or Vz 3 — z . z — z ) . (271) 

With a parameter 

v = a l + ia 3 (272) 

we should have to interchange the suffixes 3 and 0, in the first form. 

In the associated pseudo-elliptic integrals (Proc. London Math. Society, 
XXV, p. 218), corresponding to the parameters in (270), 

*x = (l — of, s 2 = c\ s s =(c-c*Y; (273; 

P = f(2 — e)(l— 2e), or f (l + e)(2 - c), (274) 

a = 2c (l — cf, or — 2c + 2c 2 , (275) 

*/— 2= 2c (l — c) 2 (2— c)(l — 2c), or 2c(l — c)(l + c)(2 — c). (276) 

Differentiating as before, and equating coefficients, will serve to determine 
the unknown coefficients H and K; the work, which is very long and laborious, 
has been carried out for me by Mr. T. I. Dewar, and he has found that 
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IP-aH* = i™^ or f-^?, (277) 

while squaring and adding in (271) gives 

K 2 + aH s =l, (278) 

whence /Tand H being determined, the other coefficients readily follow. 
Thus, for the prolate body, v = % + %ct 3 , 



E 



, _ X s — (2 — c)(l — 2c) iV^iVg 



3 



2Z 3 



so that 



^- ^+(2-^(1-20)^^3 . (2?9) 

2 fl3T=^r- (1-^(2-^(1-^^. (280) 



while for the oblate body, w = ^ <a 3 , 



m- (l +c)(2-c)N 1 N 2 -L s 
~~~ 2L 3 



^ = (l + o)(2-c)^ a -^, 81) 

2x/ 

so that 2fffir= (l + c)(2--c)-i? ^ ( 28 2) 

If we try to cancel the secular term by putting 

L = iP = i(2 — o)(l-2c), or i(l + o)(2 — o), (283) 

we find 

9-W?=(l— 2c)(l — c)(— 5 + 2c + 2c a ), or dNj = — 2(1+ c)(2— c)(l— 2c) 3 , (284) 

and these are negative for the region <C c < 1 , so that algebraical cases cannot 
be constructed. 

As a numerical verification with a parameter v = a x + ^ 6> 3 , we may take 
the case worked out in the "Applications of Elliptic Functions," p. 348, in which 

I . , _i "/72 3 — 42 + s/1 ,, , _ .„ , _. 

*-«* = i«» W 2 (!-*)■ ^(-**-2V73 + 5) 

! • _i(— z a + 2V7s— 3)* . . 

= * sin lL - 2V2(1 _ g2) l • (284a) 
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With B = and a parameter v = wj + ¥ %, 

i^ = 0, iV g = «/(l-2c), iV 3 = (l-c)V(l — c s ), 
i = V(l — 2c.l — c 3 ), if=cV(l — c 3 ) — cV(l- 2c), 
P=-f(l+c)(l — 2c), 

in which ff= 0, K=l, ^ = 0, #2 = 0, 

TT -o P - (1 + oV (l — 2c) + (1 — 2cy (1— c-j 
-^i-a n - c2 

^_ —if _ ,|V(1— o)— V(l+c.l — 2c)) 8 

With B = 0, and a parameter » = 30 3 , 

#3=0, Z = — V(l— c 3 . 2c — c 3 ), 
JV?= (1 — c) 3 (l + c), iV|=c 2 (2e — c 3 ), 
if = V(l — c 2 ) + V(2c — c 3 ) , 

in which 

#=0, jBj= — 3 ^- = (2 — c)V(l — c 2 ) + (l + cV(2c — c a ), ff 2 ~0; 

ft 

K=l, ^ = 0, iT 2 = -^^ = ^V(l-c 2 ) + V(2c-c 2 )f 2 . 

30. These algebraical calculations are too long and complicated to be 
inserted here, and the complexity compelled us, as in the corresponding work for 
the Spherical Catenary, in the Proc. London Math. Society, XXVII, to turn 
elsewhere for some clue to the values of the leading coefficients H and K, upon 
which the remainder depend by simple linear relations. 

As in the Spherical Catenary, we turn to the degenerate form assumed 
when we take 

Z = 00, u = V 3 . 
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The preceding cases enable us to infer that the general form of the solution 
for a parameter 

v = a 1 + ^, or 0!+-55- or -^ , (285) 

n n n 

where n is odd, can be written 

1 Hz™ -1 + 

X = — cos_1 — (x _ /\iV V (a — « • a — Zi), or V(g 3 — z.3 — gl), 

or \Z(a — z l .z — s 2 ) 

1 Kz n ~ x + 

= Y sin " 1 (i_^)*»* V (g» — g.g — g»), or \Z(z — g.g — g 2 ), 

or V(g 3 — g.z — 2 ), (286) 

but for parameters 



V = 6)l + _^L ) or •£»-, (287) 

n n x ' 



and n even, we must take 



1 Ez n ~ x + 

% = _ COS" 1 , 1 _ z y n 'W (z — Z.Z 3 — z), Or V(z — Z 1 .Z — Z 2 ) 

1 Kz n ~ x 4- 

= IT Sin (i — g ^W (s — 2 2 • g — gi), or V(g 3 — g • g - g«). (288) 

Squaring and adding in (286) and (288) shows that 

E 2 + K 2 = 1 , or E*—K*=(— If, (289) 

according as the body is prolate or oblate ; and putting z = <x, 
r = — tan- 1 ^^ = — sin" 1 2SfiVa 

= JL cos" 1 (5" s — JT 2 ) , or J- cos- 1 (E 2 + ^ 2 )(— 1)». (290) 

But, by means of the formula 

we may write equation (133) 

<k = » en* + * log off* + ^) + uo g ^"-^-^^7^ 

(3 (2m — «) °p (w + ©0 — p (« + «? 2 ) 

= ij>t + hil(2u,v) +klog '( u -*\-*( u - v >\ , (292) 

f> (w + vi) — 9 (u + v 2 ) v ; 
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thus adding the amplitudes and parameters of the two Elliptic Integrals of the 
Third Kind in 4i and ■in- 
putting z = cc and u — v 3 in this last relation, then from (96), (97), 

9 fa-vj-f <&-*,) f298) 

while /(2i> s , ») is the value of the pseudo-elliptic integral when u = 2v 3 or 

s ' * = p2v s — p ( Vl + »i) = w • ( 294 ) 



Jf 3 r " rv x ' " if 
s — o = L % , (295) 

8 — 8 a = aN!. (296) 

The form of i~(«) being given by 

/(«) = — sin J -^ — xnr 
v y w (« — (T)* n 

- n C0S ( a _<y)** ' (297) 

therefore, when z = oo and ^ — ^9/ is replaced by %, 

3 = iJ(20„*)+i*log(-l); (298) 

or, disregarding for a moment the ambiguities of sign, 

J_ Bin- ^y" = — sin" 1 2HKs/a, (299) 

2n L n 2n v ' 

J_ cos" 1 QJMi — _L cos" 1 (# 2 — alP) , (300) 

2rc L n 2n v " v ; 

or 2HK=^ (301) 

H % — aK*=±®*^i, (302) 

and these equations, together with (289), determine the leading coefficients 
H and iT, upon which the other coefficients depend by simple linear relations. 

Since z = o° in (167) makes & = $, we see that the leading coefficients 
land e7in (184) are determined at the same time, by taking /== K and J= H 
or /= E and J= K, according to circumstances. 

31. To make quite sure of the plus and minus signs, which are apt to be 
baffling, it is well to make a recapitulation of all the various cases which may 
occur, according as the body is prolate or oblate, and n is even or odd. 
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"When the body is prolate, and the parameter of the form 

v = % + -^ , (303) 

n 

then, when n is even, we must take 

Z = — cos * ^ (1-?)** V fa— s.fr-g) 

= — sin (1 zzjp V (g — g».g — %), t (304) 

and when w is odd, 

% = — cos (l — ?)*« — V(g»-g. g — gi) 

= ~Sm 1 (1 _#yn V (go-g.g ~g»). (305) 



But with a parameter 



where w is an odd number, 



v = o, + — , (306) 



1 .ijffg*- 1 ^- // \ 

X= — cos (1 ZZ^p V(g — g-g — gi) 

= -~ Sin " 1 %~l + dy*' ' v/(g 3 -g.g-g 2 ), (307) 

obtainable from the preceding case (306) by an interchange of z and z 3 , or 
% and z % . 

In all these cases of the prolate body 

g 3 > g > «*, (308) 

and H* + K* = \. (309) 

Putting z = oo, 

r = — tan" 1 ^ = J- sin" 1 2#ff" = -i- cos" 1 (if 2 — /f 2 ) , (310) 

and at the same time 

* = i 7(2*3, *)= ^ sin- 1 ^ = ^ CO*" 1 ^; (311) 
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sothat *HK=^r, (312) 

IP—K^— aN ff y . (313) 

Jj 



32. When the body is oblate and the parameter is of the form 
we must take 



v=^-, (314) 

n v ' 



1 _ 1 Hz n - 1 + + H n _ 1 .. N 

^=-_COS 1 (i_ g2 ^T V (%— g.a — g,) 

= i Sin_1 ^"^ ^ - TJ r^Si* - ** .* - «i) (315) 

and J2" 2 — ^= (— 1)". (316) 

Putting z = oo and it = v s , 

v = — tan- 1 § = J- sin" 1 2 J HLfiT* (— 1)" = _i- cos" 1 (# s + #»)(— 1)» (317) 

so that, allowing for the effect of the Ji log ( — l) in (298), we may put 

2HK=^, (318) 

H* + K* = ®^. (319) 

33. Finally in the case of an oblate body, when the parameter 

« = — 3 , (320) 

n v ' 

and n is an odd number, we must take 

y = J_ COS" 1 ^ + g^- 1 +....+#„ 
X W (1 — 2 2 )*' 1 

K* — H*=i; (322) 

* = 4- tan_1 # = ir cos ~ 1 (^ + ^ 2 ) • ( 323 ) 





= 


1 
n 


sin~ 


with 








and proceeding 


to 


z = 


= <» I 






£ : 


1 
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At the same time the associated pseudo-elliptic integral 

no \ 1 -1 F 1 • -i GrN a N 8 N y i /0O( x 

I(2v 3 , v) - — cos * _ = -— em ij^L , (324) 

SOthat Bz + K* = J^, (325) 

Z 

2JIff"= GN °-^N y ^ ^ 326 ^ 



In such cases it is generally simpler to take a parameter 

4g>s 



v = 

n 



(327) 



so as to make a = in (294), and now, when z = oo, 

S = Z 2 , Nl = s a -L\ (328) 

The resolution of the cubic $ in (50) is not required now, because symmetric 
functions only of the roots s a , s 2 , s 3 occur ; thus from (63), 

V— 2 



if » = ivf + m + ivi + 2i 2 — 



X 



= « x + S 2 + S 3 — Z 2 — -^ 

= i(2/+l) 2 -2a;- Z 2 - ^ , (329) 

ZJf 2 = — ay + i{(y + l) 3 — 8tc} Z - Z 3 

= JfV« + SLlPpv — L\ (330) 

as in Lame's equation of the second order. 

34. Withs = Z 2 , 

- l S={xtj + (y + l)Zy— ^(x + Uf-S^, (331) 

where S 1 = xy + 2ccZ + (y + 1) Z 2 + 2Z 3 , 

# 2 = X y — 2xL + (y + l)Z 2 — 2Z 3 . (332) 

At the same time we shall find that we can put 

i r/o \ 1 • _ x J. + 5Z + CZ 2 /e 

4 7(2» g( w) = — sin — n -1^ Vtfi 



n 



2Z* 



-^T cos 2Z^ T' (333) 
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while with z = go, 

y — — sin- 1 K= -A- cos- 1 5- , (334) 

* n n % K ' 

s „,hat Jrorg= ^±M±^+^ Vi8 , t 

thus determining the leading coefficients, when the parameter v = 4a 8 /n . 

35. For instance, with « = 3 , and taking the corresponding pseudo-elliptic 

integral, with P = \ , 

is + c 






= uos ~ ls -^= is ' m ~ 1 ^§' (336 ) 

we have K _ y(c + I? + 2Z 3 ) 

2X^ 

g.^i^ 3 ), (337) 

Then Zif 2 = c + i i - X 3 , (338) 

^b _ _ iv^, _ jv|(c+x 2 ) 2 -4i: 6 i 

il/ ZJf 2 c + II-X 3 

4c + Z — 4Z 3 ' ^ ' 

and the result is of the form 

•4/ — m* = | cos" 1 ! — ^ — » -J- — ? = I sin" 1 ^ 8 ,V ^. (340) 

1 (l — z 3 )* (l — zy 

By squaring and adding, and by differentiation, we find 

n- _ 1 + 6/ ^ _ 1 — 36Z 2 rr 

^ 2_ 2fiT Al 2JP~ ' 

3 M 2M ^ M T 

2if T if 

- 1-2X ,. yj(c + XT-4Z^ l-BL ff 
~ 2M 4c + L—4L s M 

— i *~ 2Z 4- 4 c + ^ 2 - 2X 3 1 — 6X ) K 

~ \ M f 4C + X-4I 3 M \ 

— {«-«+«(!-« >S£S=£l£- ' 3 "» 
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If we cancel the secular term by putting 

Z = iP = i, (342) 

then E 1 — , K x — , H % = o , 



and ^ x K 

' "~M 






fl" = V(54c+1), 
Z = V(54c+2), 
9if 2 =54c + 2 = J? 2 , 
3Jf= — Z", H s = — 1, 

W^^TT' ~M = 2K' (343) 

and now 

, i _!\/(54c + l)z 3 — 1 

4-, COS {l _> }i 

= $ Bin-* j—^-j^y \— (54c + 2)z« + 32»+ 2^/(54c + l) 3 — 3}. (344) 

This can also be written, with H for \/(54c + 1), 

(1 _ z y<?<H — H£ — 1 + *W| — (# 2 + 1) 2 4 + 3a 2 + 2i& — 3 \, (345) 

which by differentiation will be found to verify the relation 

^ -- 9 L — Bz /o 4fi x 

■&-~ 2 iif(i- s 2 V^' (346) 

and thus (345) represents an algebraical case of the motion of the axis OC of 
an oblate body, relatively to a line Oz fixed in direction. 

At the same time the curve of (a , (3) or (p , m) will be given by a relation of 
the form 

„ _ i cos -x & + W + jg + /. 

* - * cos (2*- ! + .#)« 



where /= Z=V(^ 2 + 1) and /= iT, (348) 

= 7z 3 + V + /,» + /, + * (Ja + «7i)\^, (349) 



or f 
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and we shall find that with 

ff-o ^-1 r- 3g T _ j_ H°+9H ,_ , . 

*'- 2 TP + l' A ~ V(H*+T)' 4 ~ 0, l3 -~(B*Tlf' Jl ~ S ' (350) 

the differential relation of (169) can be satisfied. 
The discriminant of 

Z x — — (R 2 + l)z 4 + 3z 2 + 2Hz — 3 (351) 

is — f*(fl»— l ) 3 (iP+ 1). (352) 

If iP< 1, all four roots of this quartic are imaginary, and Z is always 
negative, so that no real solution exists. 

But if IP > 1 , the discriminant is negative, and ^ has two real roots, so 
that real cases of motion can be constructed. 

36. Making use of the results of the pseudo-elliptic integrals of higher 
orders, we find for n = 5 , putting Halphen's y = x, 

jr nv rr_ (x-LW\x z +2xL + (i + x ) U+2II) 

° 2X4) or 2 ( _ L y • 

Hor K= (x ± LWW -j*L + (1 + x) I? -2U\ (863) 

2D or 2( — Z) 2 v ' 

according as L is positive or negative, and to make the solution algebraical by 
cancelling the secular term, put 

L = tP{v) = ±=£°. (354) 

This makes 

K tt— (— 1 + 13s V(3 + 33a; + 121a: 2 4- 9x 3 ) 

OT (1 — 3a;)« or'(3a; — l)« 

ffor K= (l + 7xW{2(l-Ux)(l-Ux- X ») 

(1 — 3cc) s or (3x — If ' v ' 

N X N,N % = V|( Sl - Z 3 )(s a - Z!»)(«, - £ a ) } 

_ V|2(l — 18aQ(l — llg — x 2 )(3 + 33a; + 1 21a: a + 9s 3 )} ,„... 
iooo "' l j 

jjp- 24(1 + 2x)(l - lb - x 3 ) ? (35?) 

and putting 

= J sin -. gg + g±*! + V (358) 
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Where Z= -V _!)(*_ 1 + *)- 4 (^*)', (369) 

the differential relation (46) shows at once that 

5i = and ^ = . (360) 

As the result of an algebraical verification, it will be found that with 
3a — 1 



L- — 



10 



rip- 24(2a + l)(a 2 +lla-l) 
~" 1000 

JT= 24(2a + l)(a 3 +lla— 1) 
100 (3a — 1) 

Z^_ (3a; — l) 3 

M % 24(2a + l)(a s + 11a? — 1) ' 

BL__ N t N z Ns 

IP M 3 

1_ (3a — l)W]2(18a — l)(9a 3 + 121a 8 + 33a + 3)j - 

24^6 (2a + 1) § (a 2 + 11a — 1) 

^B __ _1_ %/|2(l8x— l)(9a 3 + 121a 2 + 33a + 3)} 
M 12 (2a + l)V(a: 3 + 11a — l) 

Jj___ __ _J (3a — l) j 

M 2</6 V{(2a + l)(a 2 +lla — l)f ' 

2a 2 250a; 2 



^ AM* 3(2a + lXa s +lla— 1) ' 

fi^3 i « i ^2 _|_ 97^ 5 

(361) 



Zif 2 3(2a + l)(a 2 +lla — 1) ' 
p_ 1 ff— 6a; 3 — 181a; 2 + 27a; — 3 

and then 



3(2a + l)(a? + 11a; — 1) ' 



(362) 



Z=± 4 _ (3a; — l)(a 2 + 66a; — 1 1) 8 
S 6(2a+ l)(a 2 +lla — l) 2 

1_ (3as — l)Vj 2 (18a; — l)(9a 3 -4- 1 21a; 3 + 33a; 4- 3) \ 

3V6 (2a; + l)« (a 2 + 11a— 1) ~~ S 

1 (3a — l) 2 (26a 2 + 21a; — 21) 

18 (2a + l) 3 (a 2 4- 11a— 1) ' 

We must now take 

tt __ (13a — l)-s/(9a 3 + 121a 2 4- 33a + 3) 

~ (3a— 1)* 

K = (7^+l)^l2(l8a-l)(a 2 +lla— l)f . ( 

(3a — 1)* ' ^ ' 
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and thence we find 

TT _5_ (8g — l)\/(9a- 3 + 121a a + 33a; + 3) 

2 3 "" (2aj+l)(3a;— 1)» "" ' 
K _ J_ (12a?+ lVj 2 (18a;— l)(ce 2 + 11a; — l)j- 

3 3 (2a: + l)(3x — if ' 



(2» + l)*(3aj— 1) 

K — -A_ ^(9a: 3 H- 1 21,r 2 + 33a; +3)\/(x a + 11a— 1) 



3 3V6 (2x + l)*(3aj— 1) 
4 \/(9a: 3 H-121,r 2 +33x +3 

3\/6 (2x + l) f (3a; — 1) 

rr _10 \/(3x — l)V(9a; 3 + 121x 2 + 33a; + 3) 

4 9 ~ (2a; + l) a ~ ' 

H- X (22a; 2 + 42x + 3)Vl2(18x-l)f , . 

5 ~~ 9V6 (2a; + 1)* * { } 

A numerical verification is obtained by taking x = 2 ; this makes 



and then 



E = $, ¥=1 — E=- |; (365) 

r? 4 5 3 Vl05 5 /.->«/»\ 

Z= - 2 ~6 S - 9— S ~18 ; (366) 



and H = Wb, iT = 3\/l4 ; 

^3=5^5, K i = Wli) 
E 3 = W21, iT 3 = V30; 
ff 4 =VV5, 5^ = ^21. (367) 

37. "With w = 7 and a parameter 

« = 40b, (368) 

we take, in (49), 

x = — c(l+c) 2 , 2/ = — c(l+c); (369) 

and now we find 

lorl 

__ {c{l + c) 3 +(l+c) 2 L—I/\V\c !i (l+c) s —2c(l+cYL+(l+c—c % )L 2 +2L 3 \ 

2L i or2(—Lf "" 

HoyK 

_ (c(l+6-) 3 — (l+ C ) 2 X-X 2 }x/|c 2 (l+o) 3 +2c(l+c) 2 Z+(l+c-c 8 )X 8 -2Xn ,, 7ft v 

2Z* or 2 (— Z)* ; ^ d ' UJ 
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and the solution is made algebraical by taking 

L = $P(v) = S + 9 ° + 5( ?. (371) 

With n = 9 and a parameter 

t> = *a», (372) 

take ("Pseudo-Elliptic Integrals," p. 232) 

aj = i>»(i— jp)(i— p + p 2 ), y=p % (i—p); (373) 

and now 
Kov H= 

\p*(l—p)(l — p + p 2 ) — p*(l— p + p % )L— p % (1 — 2p)L 2 + L s \ 
*/\p i (l—p)(l—p+p 2 )+2p z (l— P )(l—p+p z )L+(l + 0+p 2 —p 3 YL+ 2L 3 } 
~ 2Z»or2(— Lf 

RorK= 

\p A ( 1 —p) ( 1 —p +p*) +p* ( 1 — p +p 2 ) L—p 2 ( 1 — 2p) X s — Z, 8 } 
^\p\l-p)(l-p+P 2 )~2pV-p){l-p+P z )L+(l+0-P z -P 3 W-2L s \ (3U) 

2L i or2(—Lf ' V ' 

and so on for the higher values of n, but the complexity increases very rapidly. 

38. Making use of the theorem in §30 that 1= Kor Rand J= H or K, 
according as the body is prolate or oblate, in the associated pseudo-elliptic 
expressions of or in (184), we can now write down the corresponding formulas, 
and determine the remaining coefficients by means of the differential relation 

d(* — pt) __ M M M y ^ ) , , 

dz ~ (1 + aE — z*)VZ ' K } 

Thus with the parameter v = Oj + •£ g> 3 , 

m — pt = i cos -1 - _^J"_ 1 s 2^( g o — z.z 3 —z) 

= * sitl_1 ! + t-z* V{z ~ * 2 • Z - Zl) > (376) 



with 



21} 



J2 = H*= L% ~ z f liV a ; (377) 
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and then 

/i ,, ^+(7+1)* 
J M 

J x _ N 1 N i -(L+lf 
T M • (378) 

The relation (376) may also be written 

{i 9 eHW ~ Pt) ) == ( Tz + 7 iM 2 o -z.z-z s ) + i(Jz + J,W(z - z, . z - Zl ) . (379) 

With a parameter « = £ a s , 

a—pt^i cos" 1 g3 _"^" + 1 ^ \/(z — Zi . Z — Z 3 ) 

= 4 8in_1 ^T^"/^ 3 - * ■ 8 - «.) . (380) 



or 

{An 

with p =H 2 = NiN* — E 



pfw-rt) J * = (/ z + ^y^ — % . 2 — z,) + i (Jz + /j )V(z s -z.z — z ), (38 1) 



22? 
J^= K l =W^^-; (382) 

and then Jj_ _ (L + P) 2 + NiN, 

J PM 

J, _ (L + P?-N,N, . . 

~T LM • (383) 

So also with a parameter 

v — % + |w 3 , or wj + lws, or io 3 , (270) 

the expressions take the form 



{£ f""-} 



3 



= (7z 3 + 7,8 + ^M^— a . 2 — a,) + * (7z ? + Jj2 + J 2 )v/(z 3 — Z.Z — Zz), (384) 
or = (/z 2 + 7 x a + 7 s )\Z(zj s — 2 . a — a,) + i (Jz*+ J x z + J a yS(zg— z . z - a,) , (38 5) 

or = (7z 2 + 7 t zj + 7 2 )\/(zj — z x . z — z z ) + i (Jz* + J x z + J z ) V(z 3 — z . z — z ) . (38 6) 
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Thus, in completion of the preceding special cases of §§28, 29, we find for 
a parameter v = 6) x + h <a 3 , and, in addition, with 

B=o,N 1 =o, L = V{l + c), 

with a corresponding 

; + V(l + c) — 1 

y(l + c )_ l 

For the parameter v=-ba 3 , with 
5 = 0, iV" 3 =0, fc»=l — E= \V(l+c)— Vc\\ I=-4 + c'), 

-, _ «* - 1 cos- 1 z - a/[> - Ml + <0-*M' l 

|V(1 + c)-VcP 

= " 8in "\^ -lvtl+ + c F^[^ [(1+2c) ^ (1+ C) " VC} "^ * (886C) 
In the numerical case of the parameter v = % + 4 « 3 in §29, we find 
n, - n* = h cos- 1 -^-8V7.+ 131 V( _ z2 _ ^ + g) 

= i sin- 1 V7g8 ( 7 3 2 l7y W7 V(— z 3 + 2V7g - 3). (386rf) 

Or, in the algebraical case, with 

£ = 0, #1=0, L — {\ — c)V(l — 2e), 

Jf= V(2c — c a ){l — c— V(l — 2c)[, AT 3 =V(1 — 2c), #3= (1 — e)V(l— c 2 ), 
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we find 

and 

11— V(l — 2c) J 
and now the result can be written 

, - pt = i COB- ^+gjffi (- * + 2 f . + t ) 

= * 3i °" $7p-$ (- * ~ 2 1 8 + ») ■ < 386e) 

JT-3J— (2-e)(l-2e ) r = _js 

1 -\it i 2 » 



where 
and then 



where 



J/ 
With a parameter v = « a +1%, and, as in (2845), 



4- sin - 1 



where r _ 3Z + (l + e)i(l — 2c) r _ ,. If 



J _ 3Z + (l+c) i (l-2c) T _ , 3 



M > ~* •• jyr _ jyr • 

So also with a parameter v = % a 3 , and B = 0, N 3 = Q, 

= * Bin " 1 ^=^iV{'-(^r^) , |' (386A) 



where #=i — .0 



_ j V(l-c a )-V(2c-c 8 ) P etc 
_ |V(l-c a )+V(2c-c 2 )j etC> 
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The expression for w in the algebraical case, when the parameter v = | o 3 , 
has been given already in (347) and (349). 

So also, corresponding to the algebraical case in (358), we shall find 

<*-*cos (p _ I + Ef 

— I a in-l ^ + ^ + Jgg + ^8 . / y 
~^ Sm («•_! + #)• VZ ' 

or 

(£^)={zr <• + #>}' 

= (/a 5 + 1^ +....+ /,) + i (Js 8 4- J^ + J^ 4- J S )JZ, (387) 



where 



/"— rr— ( 7a; + l)\/|2(18a; — l)(a; 2 + lias — 1)} 

(33 — 1)* 

r- jrr_ (13a; — 1V(9« 3 + 121a; 2 + 33a: 4- 3) ,„ ftft . 

J-H (to — 1)« ' (388) 



and then 



r _ in L y_ 5 (13a; — l)>/(9a: 3 + 1 21a; 2 + 33a; + 3) 

1 ~~ M V6 (3a; — l)V|(2x+l)(^+ 11a;- l)f ' 
_ in A r- _ J- (7^ + l)V^2(18a; — 1) [ 

I_ I V6 (3aj — l)«s/(2a; + l) ' 

T — 3 (31a; 8 4- a;— l)V{2(l8a? — 1)\ 

3_ 5 (3a; — l)V(a» + lla; — 1) ' 

7 _ l (18g— l)(31a; 2 + x — l)V(9x 3 + 121a; 2 + 33a; 4- 3) 

2 3 (2a> + l)(a: 2 + 11a; — l)(3x — 1> 3 

I ■- — 10 (122a; 4 + 309a; 3 — 3 57a; 2 + 66a; — 3) V(9» 3 + 121a; 2 4- 33a; 4- 3) 

3 3V6 (2a; + 1)»(b* + 11a; — l) s (3a; — 1) ' 
T " 1_ (31a: 2 4- a; — l)(26a: 2 + 21a;— 21) V^ 2 (18a — l)f 

3— 3V6 (3a; — l)(o* + 11a; — l)(2as + l) 1 

/ 4 = o, 

J, = -f / s (l- E) - 2JLj 3 (389) 

= ^o {d tV?Yj~™ X V) (5811W + 31956a;*- 29655a; 4 
9\/6 (2x + l)f (a; 2 + 11a; — 1)| v 

+ 37830a; 3 — 13495a; 2 + 1596a; — 57). (389) 
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The special numerical case of x = 2 makes 



I:=3\/14, /! = - 


25>/30 T 25V14 r 

- 5 , h 3 , ^3 ~ 


25\/30 
9 


</= 5/\Z5,«^ = - 


-5V21 ,/ 2 =: 35 3 y5 .4 = 


25^21 
9 


with 


if 12 ' M 12 ' 




and 


Z =: — 3* - - - Z' — -5—- 2! - 

6 9 


5 

" 18' 



r _ „ r _ 25V30 
4-0, i 5 jg~ 



and these numbers will be found to verify. 

39. In some cases a factor of Z, say z — z , can be obtained; and with 

Z=aa 4 + 4&z 3 + 6cz 3 -f 4cfe + e, (27) 

we can take "Weierstrass's formula 

pu = i (oag + 2fe + c) + ggo±M±J£gl±J , (390) 

z — z 

and now we put 

—jj^r — V ~P U ' ( 391 ) 

where <pv has the expression given in (86). 

The more general form of T, involving circulation, etc., may now be 
employed, without introducing additional complication. 

It is convenient to put 

l + aE = &, (392) 

so that 

Z = a (z* - l)(rf - V) - 4 (^^)"; (393) 

thus a =F 1 is a factor of iT if 5 = ± Z, and z =f & is a factor if B = ± Z&. 

The case of B — L may be supposed to be produced when the body, pro- 
jected originally as a perfectly centred projectile from a rifled gun, has struck 
an obstacle, thus setting up gyrations in which the axis OG periodically passes 
through Oz, as in the corresponding rosette curves described by the axis of a top 
discussed in Professor Klein's paper on " The Stability of a Sleeping Top," in the 
Bulletin of the American Mathematical Society, Jan. 1897. 
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40. Now, with B = L, and z — la factor of Z in equation (390), 

pu = i (a + 6) + " + *» + g , (394) 



c = -|«-|J-l4, (395) 



and ! ! ro „ i? 

If 8 

rf=2 § =2 y ; (396) 

«othat p „_ pM = i ^ + _L» , (397) 

z — 1 

a — s = iE Z+ * ; (398) 



or from (65), 



M 2 4 2 — 1 ' 

cr — s a 1+2 

V— 2 "2X 1 — z' 

c — s a _ _a_ 1 + g a 

V— 2 2Z 1 — g.' 

* — s„ . _ a 3 — z„ 



V-S i (1 -*„)(! -a) ' 



But from (393), 



if 3 (1 - Za)(l - *) ' 



so that 



and i S _ a m Z 



■E 



M* -.a(l-z a )(l-z,)(l-z y )(z-iy 

if - ¥ ■* ( a _ i)* ' 



(399) 
(400) 
(401) 



° r S ~ S "= IE-. ^^ r. (402) 



a (z — z a )(z — z,)(z — z y )=a(z + l)(a a — ¥) — 4 ~ (z — 1) , (403) 
a (1 — 3 a )(l — Z<J )(1 — as,) = 2a (1 — *») = — 2E, (404) 



(405) 



and ds __ t „ dz , 

w-* E jjn=r?' (406) 

so that, as in (48), 
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With L = B in equation (59), 

J L ^^B_^_ a N 1 N i N 3 
M M LM* ' 

or DM = aN x N % N z , (408 ) 

so that, from (63), 

ai 4 iP = i 4 (l? + SlPfv - a ^J-^) 

= aN?NiNi = (I? + <s — Sl )(L* + a — s,)(I} + a - s 3 ) 

= Z 6 + SL 4 M 2 pv + j- i 2 ifV'« + i 2 , (409) 

or aiV— 2 +%I}MYv +12 = 0, (410) 

so that L must be determined from this cubic equation in a very similar manner 
to that required in the Spherical Pendulum, as discussed in the Proceedings of 
the London Math. Society, vol. XXVII, p. 607. 

Having determined L, or, in a pseudo-elliptic case, having expressed 
L, a, a/ — 2, pv, p''v and P in terms of a single parameter, 



^ = 2 Z, l dz 



M 1 + s s/Z 
2L ds 



(411) 



so that 



1 + z V/S" 
and 2 + 1 a — s L , x 

a^i ~ iEM 2 ~ 1«V- 2 {a S) ' 

z + l = T7Z 2L l a ~ S \ , t (412) 

i (c — s) — i a\/ — 2 

, i Z(<r — s) — i«\/ — 2 <fo 

= (i-iP)^ + i^^l^V=2 A (413) 

4> = ^=^nt + /(„), (414) 

and similarly for t», which can be made to depend upon the same integral I(v). 

The case of B = — X is the same as the preceding, with z changed 
into — z, so that this case does not require separate treatment. 

A reference to equation (14) shows that, with B = ± L, 

#1 = or *b = 
dt U ' ° dt U ' 
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so that d<p /., C\ , diL / A -,*\ 

and ty is now pseudo-elliptic with •fy, as also u + vi and £> + qi, on reference to 
(179) and (180). 

41. A numerical case will serve to elucidate the preceding theory. 
Take v = f o 3 , and in (410), a = — 1„ ?/ = cc = — J, when 

*/— 2 = a: 2 = i, MYv = x* — x s ~%, (416) 

and then — i £ 3 + T V Z 2 — A = , (417) 

which is satisfied by L = — i; and now, in (353), 

H=2aJ2, K=—3. (418) 
Also in (329), 

M' z = i{l + x) 2 — 2x — L z — -S-=2, (419) 

ZJf 2 = -i, ~ = -|V2; (420) 



and in (65) and (59), 



BLM= - A^M = - 3^2 , (422) 



1 z 



*T+V 



and now in (147), 



(423) 



Jtf ~~ 4 V 2 ~~ M ' 

and Z = — z*(z* — 1)— 1(3+ l) 2 

= (z + 1)(— z 3 + s 2 - iz - i) . (424) 

We therefore take, from (398) and (405), 



(425) 



« + * = r+v ( 426 ) 

^(ifzF <« 7 > 



-£- = — iV2; (428) 

n 
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so that 

^ =iV2 \ (429) 

dz (1 — z)\/Z 

d^-pt) ^ -z 80) 

dz (1 — z)s/ Z 

Taking the pseudo-elliptic integral (Proc. L. M. S. XXV, p. 214) 



T( i ,. \ _ i / ^(l — 3z)g — a 



vs ds 

_ x t (1 - 3a;) a»— (2s 3 - x s ) 8 + * 4 
- T cos — t 

= i S in- 1 ^^V/S', (431) 

with y = a- = — i, and making use of (425), we find that 

, . . _! 2V2z ? — |V22 + f-v/2 ,, , ,. 

^fsin-^ -^ + V C-^ + ^-^-l), (432) 

(1 — zf 

and this, on differentiation, will be found to satisfy (430). 
As for equations (165) and (169), they become 

It 

An 



(433) 



and ^=W2iil, (434) 

dz — zs/Z 

0r gfeL=lg) = iva *, ,; (435) 

the integral of which, derived from the pseudo-elliptic integral (Proc. L. M. S. 
XXV, p. 213), 

772 „ \ _ i ^* l(3 + g)(a -j- as) — x 

.! (3 4- a;) s 2 — (1 — 4a — 2a; a )s + a; 8 + x 3 

1 2 (s H- a;)* 

= isin- l g ~^ + fV ^, (436) 

2 (s + x) 5 



= | cos" 
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by making use of (426), becomes 

« — pt = | cos- 1 ? s JfA*^ii V(z + 1) 

= I sin" 1 ^±ilv(- 2 s + ^ — iz-J), (437) 

(— s ) 5 

which will be found to satisfy (435) . 

Equations (432) and (437) can also be written 

]tan 8 i^-P" i j«=:V2 ^.""j*. t f + * ~ 3 f t V~ g " "VjT^ Z ~* ' ( 438 ) 

j^pe^-^ 4 } 8 = (a a -|z— i)V(z + l) + *V2(z + |V(— s s +2 3 — iz — *). (439) 

Now, in (179), 

P(u + vi) =—Fsm e -( 1 -!) rt * + *" i e «'-^x, (440) 

while 

{sin e^-^'p = \/2 (2z 2 — 1 z + |)(s + l) 3 + i (— 3s + l)(z + l)VZ. (441) 
But, while keeping Or unchanged, G and r may be varied so as to make 

(i — x) r = ^' (442) 

and now 

j P (u + m)l *= — V2 (2a 2 — ? z + |)(z + l) 3 + * (3z — l)(z + l 2 )\/^> (443) 

with a similar expression for^> + qi obtained from (180). 

42. With B — Lh = Q (444) 

and l + aE=*r, (445) 

'=(-t)H'-<' + x>- 4 £(-4)1- (446) 

A similar calculation shows that 

- 1F - = ia(- T r-l) F - («7) 

a- x 



and aZV— 2 — 1 £V"fl — i 2 = , (448) 

the same as (410), when L is changed into — L. 
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/? d R 

With z jr- a factor of Z, — j— = when z = - T - , and the curve (a , /3) 

Jb dZ h 

has a series of cusps. 

But with z 55- a factor of Z, -^- = when z = -=- , and the cone des- 

H dz M 

cribed by 00 round 00 has a series of cuspidal edges ; and in this case we find 

ti> u—pw = ia(~—l') r?.. (449) 

Z ~ ~lT 

43. The integrals are non-elliptic when the quartic Z has a pair of equal 
roots ; this will be the case if the body is projected in the direction of its axis 
OG, with a rotation about 00, as if fired from a rifled gun and perfectly 
centred. 

Denoting the velocity of projection by V, and the angular velocity by r, 
then in the preceding notation, 
F=BV; G= GtRV= OrF, or B — L; E—0; 

2T= B V" + Or 1 - Jg- + Or* ; (450) 

Z=a^-lf-4^(z-lf 

= (z-iy[a{z + lY-4^, (451) 

ndt 1 

^ ~ m „wJ„„. „v TTfl ( 452 ) 



(i- a y|«(i + 2 ) 2 -4^J I 



dz M 



(-^{„, 1 + „'-4|| <« 3 » 



If the body is oblate and a = — 1 , the only solution is a = 1 , so that this motion 
is stable ; but, at the same time, putting 



1— z=tf, (454) 

_ — *,»»* ( ' . 
df 



then ^- = «V (- 1 - ^ + tf - i y 1 ) , (455) 
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so that ^f-_ B «f 1+ & \ v , 

~-« 3 (l + ^)2/. (456) 

when y is small ; so that the number of complete oscillations in one second is 

£V0 + £) • <« 7 > 

But in a prolate body, with a = + 1 , then (i) with -=^ > 1 , 

dt 1 



^ rt _jihx. w _-pr ( 458 ) 



and integrating, 



= ^ 3^-, /(jr- 1 X 1 + - | - , y) (4M) 

//i 2 \ / r ' ( 45y J 

V(^-l) V 2^(1-.) 

but this does not represent any real motion of the axis, since z oscillates 

between ± 1. 

The only solution is therefore z = 1 ; and to find the period of a small 

oscillation, putting 

l — z — if (454) 

as before, 

^ = K y(l-^-y + ^«), (455) 

v=- n '(v- i )y + < 466 > 

and the axis thus makes 
complete oscillations per second. 
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Thus to secure the stability of this elongated projectile in its flight we 

must make 

L Cr 



M 2An 
44. But (ii) with -=-= <C 1 , integrating 



>1. (458) 



nt = —, ^5— ch 






= ! sh-i / ± m '^ ILL, (459) 

so that the motion of the axis is unstable, and, after an infinite time, is given by 
the equation (459) above. 

Putting L _ n ^ (ARn\ 

° -^=cosa, (460) 

i «/ . . \ (1 — cosa)(l 4- z 4 2 cos a) /ahi\ 

ch* (nt am a) = ± — ^— — —. -, (461) 

v ' 2 cos a (1 — z) 

,», , . x (l+cosa)(l+s — 2 cos a) (aro\ 

sh 2 (nt sin a) = > — ! - — ^ — j- . , (462) 

' 2 cos a (1 — z) 

u fn j. • \ z — cos 2a /,<<jq\ 

ch(2«i5sina) = t-—. — r, (463) 

v ' cosa(l + s) 



cos 2a + cos a ch ( 2nt sin a) 
1 -f cos a ch (2nt sin a) 



(464) 



i i o cos a + ch (2nt sin a) (ark\ 

14a=2cosa— -J- , v - . (; (465) 

1 4" cos ach (2n£ sin a) 

d^i 2n cos a 

1ft 1 4-s 

1 4- cos a ch (2nt sin a) 

~~ * cos a 4" ch (2nt sin a) 

i wsin 2 a [Aaa\ 

= ncosa-\ , t- la . • v , (466) 

cosa + ch(2n^ sin a) 

. . , , _i 1 4- cosach(2ntf sin a) 

4- = nt cos a + i cos ; — ,/ OJ . e*- 

cos a 4- ch ( 2re£ sm a) 

. ,i_i 2 cos a tAan\ 

= wi cos a + J cos 1 — — . (467) 

1 4- z 



Greenhill : Motion of a Solid in Liquid. 63 



At the same time, from (169), 



dm L z — i 1 _ z 

■ = 2n -Tj-r 5 = 2w cos a 



dt M 1—Z Z ~ 1 + 2 

cos 2a + cos a ch ( 2nt sin a) 

— ^2, - - . 

cos a + ch ( 2nt sin a) 

w sin 2 a . . 

= ncosa xmm—- r > (468) 

nrvs rt on I vmr sin « l v i 



cosach(2wtf sin a) 

„ . i „„ _i 1 + cos ach (2nt sin a) 

™ = nt cos a — cos * — ! — = — i : J- 

cos a + ch (2nt sin a) 

, 1 _i 2 cos a 

= nt cos a — i cos 1 ; 

1 +s ' 



(469) 



(470) 



so that, as in (167), 

^-^cos-i^^^sin-^l + ^-^osM 

1 + Z 1 + 2 

Since JE= 0, a= 1, 

so that the (a, /3) or (p, m) curve is given by 

1 _i 2 cos a 

m = ra cos a — -J cos 1 ■ , „ 9 ., 

= w< cos a — 1 cos -1 ^ — -2 ■ /±ay > m (47 2 ) 

J. 3 W 2 

Prom equation (155), 

_ ^ , jjj, 2 sin 3 q ) 2 

■" P "•" t l+cosach(2^sina)r (473) 

or efy/ y . An 2 sin 3 a 

e& .P 1 + cos a ch (2w£ sin a) 

, . An? sin 4 a 

"*" F { 1 + cos a ch (2n* sin a) } ' (474) 

and integrating, 

y An sin 2ash( 2^ sin a) 

r ~~ ~F 1 + cos a ch (2^ sin a) ' ( 475 ) 
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thus determining a, /3, y as functions of the time t, an infinite time after the 
start. 

If Z = 0, cosa = 0, a = in, (476) 

as if the body was fired from a smooth-bore gun without rotation ; then 

d^ _ n dra _ 

while efe ,, . „,. /^.7-7\ 

_ = »(l + a»), (477) 

a=tanhw£, (478) 

p = -^ r sechn^, (479) 

and the motion is confined to one plane. 

45. When the body is projected sideways, in the direction OA, with velocity 

V, and with component angular velocities p and r about OA and 00, then 

initially, 

u=V, v — 0, w=0, 3=0, (480) 

so that F=PV, G = PVAp, 

^-=A^ = 2 4; (481) 

n AnF M ' K ' 

2T—Cr* — ^ = PY % + Ap z + Or 2 — *L. J % — Or* 
Jti K 

= -F*(^-^)V* + Ap> 

= — An 2 a + 4 : An i ~; (482) 

J9 =-« + 4 i^' ^=-« + 4 ^; ( 483 ) 

= B j«(^-l)-4^^ 8 + g^} ; (484) 

so that this case comes under the head of those cases where a factor of Z is 
known ; and here we have to put 

a — s 



M* 



■pv — pu 
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When there is no rotation about 00, so that r = 0, jE?=0, 

Z = 3*{a(**-1)-4^|, (486) 



ndt 

dz 

z 



V / {a(^-l)~4^|, (487) 



and therefore z = is the only solution for a prolate body, with a = 1 , since 
a 2 — 1 is negative ; but 

so that the axis makes 

^A i + *£) < 489 > 

small oscillations a second. 

But with an oblate body, and a = — 1 , 

eft 1 



dz .If, ,V ■ .V (490) 



V0-«£-') 



i 3 

so that the form of the integral is different according as 1 — 4 -==- is positive or 

negative. 

When 1 — 4-^ is negative, the only real solution is again 3=0, with a 
small nutation of the axis 00, making 



n 

oscillations a second. 

L % 
But, with 1 — 4 ^p- positive, 



•l^w*- 1 ) < 491 > 






tttf = _^ -__ c h 



Z = 77 4^~- (492) 
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Also 

<ty _ 2 L 1 



dz M ,, ^. //, , L 






M dz ^ M 



(I-VCI-^-*)' (493) 



2 A 

i=24«<+ sin- 1 i¥ 



M \/(l — s 2 ) 



.A*-**-*) 



With 29=0, 27=1, a = — 1, and in (165), 

|£ = *-l+tf=*, (495) 

p = ^ , = ^/(l - 4 £) S ech sj(l - 4 £) ««, (496) 

and from (169), 

4^ = 2A l- 2 LruM (497) 

<fe Jf \/^ M dz' K ' 

ns=2^ra, (498) 



i!f 
so that 



An jf, . Z a \ , y\}~ A W) 



a curve of the nature of a separating herpolhode. 
At the same time, in (155), 

dt P 



M 



(499) 



__7-^(l-.4^) S echV(l-4^-) B <, (500) 

3 



eZy pr _ An 

~W~ ' ~F~ 

r = Vt-p^l-4^) tanh ^/(l - 4^)nt, (501) 

thus determining the motion of the body completely 
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46. The discussion proceeds in the same manner when we put 

p = 0, G = 0, Z = 0, (502) 



and now 



Z—z 



dnt 
dz 



'{<*(**- 1)- * |lj, (503) 

V{«(,-.)-4£} ' ^ 



d *= — 2 B 



da M (1 — z 2 )\/Z 

B l 



= — 2 



m ,, -./L^ tx A &V ( 505 ) 



(l_^ja(«i»-l)-4^.| 



and since 



D2 

#=-«+• 4 ;pr> ( 506 ) 



dz M (1 +aE—s?)*/Z 

B 



= 2 



i¥ /. £ 3 A / \ .. ^ . B* I • (507) 

Considering only the case of the oblate body, a = — 1 , with 1 — 4 -=p- 
positivo, 

Z = ^(l-4^) 8 ech^(l-4-^-) n (, (508) 

■4> = cos * 



A 1 - " ifM 1 - 2S ) 



— oin - 1 - 



= sin - ' -77 wr (509) 



m = cos 



-1 



z 



= sin" 






(510) 



34 = 4 y +r; (511) 
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so that 



</('-«#) 2— -M i -3). <•»> 



... -4n 
or, putting ~y = o, 



a 



M 2 (3* __ 



g + * ^2~ "^a - = 1 > ( 5 1 3 ) 

the equation of an ellipse. 

The expression for Fwill be of the same nature as in (475). 

This last state of motion can be illustrated experimentally by holding a disc 
of cardboard between the finger and thumb and flicking it sideways in the air. 

47. The analysis developed here for the motion of Solid of a Revolution 
under No Forces in Infinite Frictionless Liquid is also available for other 
dynamical problems ; for instance, in the exact treatment of the Precession and 
Nutation of the Earth's Axis, when, as in Poinsot's "Addition a la connaisance 
des temps," 1858, the mass of the disturbing Moon and Sun is supposed dis- 
tributed in the form of a circular band or ring in the plane of the ecliptic, or 
else condensed into two equal repelling spheres, placed at the opposite poles of 
ecliptic at the same radial distance. 

These equations are considered by Tisserand in the Comptes rendus, t. 101, 
1885 ; as also Gylden's intermediate orbit, described under a central force vary- 
ing partly as the distance and partly inversely as the square of the distance, 
which leads to similar analysis. 

The motion of a ball rolling on a gravitating sphere, in which variations of 
internal density give rise to zonal harmonics of the first and second order in the 
expression of the potential, as well as the motion of a particle sliding on the 
smooth surface of a homogeneous gravitating ellipsoid of revolution, are further 
applications of the same analysis. 

48. Taking the case of Precession and Nutation, in which the disturbance 
is due to a mass M, distributed in the form of a ring of radius B, the forces 
acting upon the Earth, an oblate spheroid of which the equatorial and polar 
moments of inertia are denoted by A and G, are equivalent to a couple round 
the line of nodes, of moment 

*ZJ£-{C— 4)sin$cos£, (514) 
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tending to decrease 3, where $ denotes the obliquity of the ecliptic, and y the 
constant of gravitation ; so that 4* denoting the longitude of the node, the equa- 
tions of Energy and Momentum may be written 

iA d ^ + iA S m^^ = -i^.(G-A) S m^ + ff, (515) 

_4sin a £-^ + Or cos §=G, (516) 

where H and K are the constants of the problem, and r the constant angular 
velocity of the Earth about its axis. 

Putting cos $ = z and eliminating -~ between (515) and (516), 
£ = -,£ ^ (l V + ¥(W)-(^) 2 

= ^|_ ( ^_ 1)(a2 _ 1 + z) ) -(^^) a | 

— n % Z, (517) 



on putting 



and then 



*&£^A = * ^=n»D, (518) 



B 3 A 'A 

d4> _ G—Crz 
dt A sin 2 $ ' 
d4f 6r — Grz 



(519) 

~dz"~ An(l—*)VZ' (520) 

equations of exactly the same form as those required for the motion of an oblate 
solid in liquid, so that the previous pseudo-elliptic solutions are immediately avail- 
able; for instance, the algebraical solutions given by equations (345) and (358). 

To gain some idea of the actual magnitude of the constants in this problem 
in the actual case of the Earth, let p denote the mean angular velocity of Pre- 
cession and a the mean obliquity of the ecliptic, then (Quarterly Journal of 
Mathematics, vol. XIV, p. 173) 

_ 3yif C — A . . 

^~ W ~~ tfr~ C0S(a ' < 521 ) 

n* — 2 -j- r(i sec a ; (522) 
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and in this equation we may take 

-j- = l, o=23°27', (523) 

and, with a year as unit of time, and an annual Precession of 50". 25, 

K(\ OK 

r=2*X366, 1* = ^*; (524) 

thismakes „= 1.115 and ^ = 2062. (525) 

But if we try to utilize the algebraical case given by equation (345) we find 

from (343), 

Gr _ B _ H _ H . 

which, with i?">-l, must be < 0.707, so that the angular velocity of the Earth 
would require to be reduced to about one three-thousandth of its present amount 
for this motion to be possible ; and now the path of the pole in the sky would 
give rise to interesting astronomical speculations. 

49. In Gylden's orbit, with polar coordinates p =— and m, the central 
force 

p =w + -4-'> ( 527 ) 

or, more generally, 

P=ap + 2S — -M-, (528) 

r 
suppose ; so that 

and a dm -, ,~ nn \ 

P^=h. (530) 

Thence 

p 2 ^- = ap* + 4&p s + Sep 2 + 4dp — h\ 

= B, (531) 
suppose; so that 

dm h dm — hu , . 

dp~ ~~ p^/R ' ° r du ™ V(— W+4c^ + 6c« a +4&w+ a) ^ ' 
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is the differential relation for the orbit, and thus pseudo-elliptic cases can be 
constructed, immediately from Abel's analysis (GEuvres completes). 

50. For the motion of a ball, of radius b, rolling on a sphere of radius 
R — b, under a zonal potential denoted by V, the equations of Energy and 
Momentum become 

0+ £)■""# +&"»•»=*• (633) 

where h denotes the radius of gyration of the ball. 

The variable part of V depending on cos 3 and cos 8 ^ or sin 3 ^, we can 
bring our equations into the same shape as before by putting 

V+ H=ln*I&(l + ~\a sin 3 $ + 4b' cos^ + D), (534) 

<®£ + S i n 2 3 ^ = n z (a sin 3 3 + 4V cos 3 + D) , (536) 

Bin»&-^=*». y • (537) 

Eliminating — ^ and putting cos $ = z, we obtain as before 

dz _ . r? 
df~ nVZ ' 

where Z=a(l-z*T + (4b>z+D)(l-z»)-4 (^=g^)' (538) 

and <H> _ o L — Bz . . 

lk- 2 M(l-z*WZ' (539) 

of which the special case for V = has received consideration, and the solution 
proceeds as formerly. 

When the ball is projected without any spin, so that r= 0, the equations of 
motion are the same as for a particle on a smooth sphere, in the same field of 
force ; and now, with b' = , 

Z= a (1 -«■)(!-«? +aD)-4^g. (540) 



so that 
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With a positive and equal to + 1 , we can put 

dz> 



df -n'Ca'-^CS"-^), (641) 

and taking 

/?>l>a>2> — a >-l>£, (542) 

the particle or ball crosses the equator, and 

nt= r <k 

=T m "(^--fy <643) 

or z — a sn fint. (544) 

Putting z 9 = 1 makes 

4^ = (/? 2 -l)(l-a 2 ), (545) 



(546) 



so that d$ __ Vjffl-l)(l-a 8 )| 

<fe (1 — 2 2 )V{(a 3 — *){fP — s 2 )}' 
Now if we put 

Ar = x andA- = a, ^=^-5 (547) 

, _ i / V{ — a(l — a)(l — go)] & 

V a — a? " """ Vja;(l— ,*)(1— /£ 2 x)} 

where X=x(l — as)(l— Fee), (549) 

and A is the value of X when a replaces a; ; this is a canonical form of the 
Elliptic Integral of the Third Kind, with the appropriate multiplier V( — A) in 
the numerator. 

With a negative and equal to — 1 , we can have 

Z=(a 2 — z 2 )(^^/3 2 ). (550) 

With Z = (a 2 — s 2 )(^ — /? 2 ) , (551) 

the path of the particle is bounded by two parallels of latitude on the same side 
of the equator, and 

l>a>3>/3>0> — /?>«>— a>l; (552) 



and 
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nt -C <& 

2=adn ant, (554) 

(555) 
(556) 



4 — - 

M* 


= (1- 


-a 2 )(l- 


-/3 1 ), 




d± 




V{(1- 


a 2 )(l — 


•^)| 


dz 


(1- 


-aVU** — ^0 


^~W 




2 S 

a? 


1 

■ 3/, —s- — 

a 


= a, 1 - 


-J 3 '-*' 

2 K ' 

a 



P^ting „ _„ . _„ , ^ __„ (55?) 

makes ^=4 / m— 4 M* r558) 

where JT= je(l — sc)(a; — &"). (559) 

With Z— (a? — z*){z % + (3 2 ), (560) 
the particle again crosses the equator, and 

, _ P dz 

1 f z a, ) 

= V(a 2 + ^) Cn " 1 I IT* ^V+F)f' (561) 

g = a en V(a s + /3 s ) «<, (56 2) 

4^J- = (l-a*)(l+^), (563) 

and with z % 1 a 2 72 /v«^ 

^ = aj - ^ =a ' XT?-*' (564 > 

^ = *f / ( -f ) / ^ , (565) 
t/ (a — £c)V-£ 

X=x(l — x)(x+ -Tj-J. (566) 



where 



Thence pseudo-elliptic cases can easily be constructed, but the secular term 
cannot be cancelled. 

10 
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But when the ball is projected with appropriate spin r, we have an addi- 
tional constant at our disposal, and it is possible to make the path an algebraical 
curve ; thus, for instance, from (345), the equations 

sin 3 0cos3^=#cos 3 — 1, (567) 

sin 3 0sin3^=cos0Vj — (# 2 + 1) cos 4 + 3 cos 3 + 2i7cos0— 3}, (568) 

represent a possible algebraical trajectory of the centre of the ball, rolling on 
this gravitating sphere, the circumstances of the initial projection being chosen 
appropriately. So also equations (358). 

51. In the motion of a particle on a smooth, homogeneous, gravitating ellip- 
soid, bounded by the surface of revolution about Oy of an ellipse 

^ + Jf =1) ( 569 > 

referred to Op and Oy as coordinate axes, the variable part of the potential on 
the surface may be put equal to \ Ap 2 , so that, with cylindrical coordinates 
p, to, y, the equations of Momentum and Energy may be written 

t'^=X, (570) 

From (569) we find 

dy (3 p dp 

"dT~~oT V(a 2 — p 2 ) dt ' 
so that 

df p 2 ^ a? a 2 — p 2 df~ P ^ ' 



(572) 



or 



\a*-(l-§-)p*\p 2( ^ = (a 2 -p 2 )(Ap* + Bp*-K 2 ). (573) 

Putting p s = a?z, 






Fo%f (674) 



z 
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where 

*=(.-l)j(l-£) — lf(^+£.-£) (575) 

'£ = #• <-) 

so that / /? a \ 

^ = 4-^ \ "?/* (577) 

tfs a 2 «VZ 

and this can be treated in the same manner as the previous equations (531) for 
Gylden's orbit. 

Artillery College, Woolwich, 27 April, 1897. 



